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Any physical process can be represented as a quantum channel mapping an initial state into a final 
state. Hence it can be characterized from the point of view of communication theory. Recently a lot 
of efforts have been devoted to encompass memory effects within this approach. It is usually meant 
that memory effects arise when the application of a channel on several inputs is not merely the 
result of identical and independent maps. Alternatively, if the channel is defined by its infinitesimal 
generator in a continuous-time description, one refers as well to memory effects if such an evolution 
is non-Markovian. The consideration of such effects has given rise to more powerful tools than the 
traditional quantum channel approaches for describing quantum information processes. Moreover 
it constitutes a fertile ground where interesting novel phenomena emerge at the intersection of 
quantum information theory with other branches of physics. Recent developments in the field of 
quantum channels and memory effects are here reviewed. 
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Any physical process implies a state change. Hence the 
"channel" formalism, commonly used in communication 
theory to map input to output information, is suitable 
for describing also physical processes, provided the input 
(initial) and output (final) states are regarded in terms 
of their information content. As consequence physical 
processes, viewed as channels, can be characterized by 
their information transmission ability. The seminal work 



of Shannon (Shannon 1948 Verdii 1998) showed that a 



memoryless noisy channel, where the noise acts indepen- 
dently on each symbol sent through the channel, can be 
parameterized by a single quantity, the capacity of the 
channel. Shannon defined the capacity as the maximum 
rate at which information can be reliably sent through a 
channel, with the use of suitable coding, and proved its 
expression in terms of entropic functions of the input and 
output systems. The work of Shannon has been extended 
to include channels with "memory" , where the noise is 
no longer independent and identical on different inputs 
(IGallagerlpesl). 



Nevertheless, quantum information theory has taught 
that it is the laws of physics that establish the ulti- 
mate limits in processing and communicating informa- 
tion (iNielsen and Chuangl 120001 iPetzl 120081). That is 



why quantum mechanics, the most advanced physical 
theory, comes into play. For quantum channels, several 
notions of capacities have been introduced, depending 
on whether classical or quantum information is conveyed 
and if additional resources, such as pre-shared entangle- 
ment, are present. A straightforward extension of the 
Shannon results leads to the so-called entanglement as- 



sisted classical capacity ( Adami and Cerf 



1997 Bennett 



et al. [1999 ), i.e., the capacity in transmitting classical 



symbols encoded into quantum states assisted by pre- 
shared entanglement between the sender and the receiver 
( Bennett and Shor 1998 1 . The classical capacity of quan- 



tum channels without any side resource has been proved 
to be expressed in terms of the accessible classical infor- 



mation (Holevo 19981 Schumacher and Westmoreland 



1997). Additionally the quantum capacity, i.e., the ca- 



pacity in transmitting quantum resources as coherence 



and entanglement, has been considered (Devetak 2005 



Lloyd 



Whi 



19971 



. e the majority of the work in quantum chan- 
nels has been initially concerned with the study of maps 
that do not exhibit memory effects, the assumption of 
vanishing correlations — or memory — in the system- 
environment interaction is not always justified. For ex- 
ample, with increasing transmission rates in communica- 
tion channels, successive transmissions happen so rapidly 
that the environment may retain a memory between in- 
dividual transmissions. Optical fibers are an example in 



which such effects may be explored experimentally ( Ball, 
Dragan and Banaszek 2004 Banaszek et al. 2004). In 



quantum information processors, especially in solid state 
implementations, qubits may be so closely spaced that 
individual elements of the environment will interact with 
several qubits, thus leading to cross-talks and correla- 
tions in the noise (Hu 2007). The consideration of spatial 



and temporal memory effects is therefore becoming in- 
creasingly pressing with the continuing miniaturization of 
information processing devices and with increasing com- 
munication rates through channels. 

In recent years an increasing attention has been de- 
voted to quantum channels with memory. Apparently, 
the interest towards information transmission through 
quantum channels with memory spread after the model 



introduced by ( Macchiavello and Palma 2002). Mac- 



|chiavello and Palma 2002 gave an example of a qubit 
channel with Markovian correlated noise in which the 
encoding of information by means of entangled input 
states may increase the transmission rate of classical in- 
formation. Such an increment in the communication rate 
due to entanglement is strictly related to the property of 
superadditivity of quantum channels for the transmis- 

Subse- 



2009) 



sion of classical information (Hastings 
quently, the study of quantum channels with memory has 
largely been confined to channels with Markovian corre- 
lated noise with the aim of deriving bounds on the classi- 



cal capacity (see e.g. (|Bowen and Mancini| 


2004 


Bowen, 


Devetak and Mancini |2005| Hamada| 2002 


). Then, cod- 



ing theorems have been devised for wide classes of quan- 
tum memory channels having structural properties that 
guarantee "regular" asymptotic behavior (IBjclakovic and 



Bochel 



Werner 



iorjBji 

2008, Datta and Dorlasl [20071 [Kretschmann and 



2005). This fact can be traced back to the work 



of (Dobrushin 1963), who showed that a limiting expres- 
sion of maximal normalized mutual information holds for 
a wide class of channels exhibiting stationary or ergodic 
behavior. A completely different approach was taken by 



Hayashi and Nagaoka (Hayashi and Nagaoka 2003), that 
applied the "information-spectrum" method to obtain a 
coding theorem for the classical capacity, following the 
work by [Verdu and Han 1994 on classical channels with 
memory. 

A lot of work has been undertaken also along the di- 
rection of bosonic Gaussian memory channels. Inspired 



by the seminal paper of Shannon (Shannon 1949) on 
classical memory channels generalizing the formula for 
the capacity of (power-constrained) white Gaussian chan- 
nels to dispersive/nonwhite Gaussian channels by the 
"water-filling" formula, quantum "water-filling" formulae 
have been recently derived ([Pilyavets, Lupo and Mancini" 



2009 Schafer, Karpov and Cerf 2009). Moreover, rig- 



orous justification of the water-filling formula for clas- 
sical dispersive/nonwhite Gaussian channel s usually ap - 
peals to the Toeplitz distribution the orem (Gray 1972 1. 
Then, the latter has been used in (Lupo, Giovannetti 



and Mancini 20101 for evaluating capacities of a general 
bosonic Gaussian memory channels that can be "unrav- 
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elled" into memoryless ones by normal mode decomposi- 
tion. 

Aside from that, models of memory channels were in- 
troduced in which memory effects arise from the interac- 
tion with an environment in a correlated state ( [Giovan 
netti and Mancini 2005), allowing one remarkable links 



between capacity of quantum channel and many-body 
physical properties of the environment, like phase tran- 



sitions (Plenio and Virmani 2007 2008) 



Moving towards the dynamics of open quantum sys- 
tems it is worth noting a Lindbladian approach to mem- 



ory channels taken in (Daffer, Wodkiewicz and Mclver 



2003[ [Daffer et al.\ |2004f ^ Then in the continuous- 
time description of open quantum systems consider- 
able attention has been devoted to the study of quan- 
tum channels arising from non-Markovian dynamics. A 
different terminology is commonly used in this frame- 
work, since Markovian-correlated dynamics correspond 
to the memoryless setting, and dynamics with mem- 
ory are associated with non-Markovian effects. A typ- 
ical approach to the dynamics of open quantum systems 
uses the Nakajima-Zwanzig projection operator tech- 



nique ( 


Breuer and Petruccione| 


2002 


Nakajima 


1958 


Zwanzig 


1960 1 which shows that under fairly general 



conditions, the master equation for the reduced density 
operator takes the form of a nonlocal equation in which 
memory effects are taken into account through the in- 
troduction of a memory kernel. Then, the problem is 
to find those conditions on the memory kernel ensuring 
that the time evolution map is a bona fide quantum chan- 
nel (Chruscinski and Kossakowski 2012). Contrarily to 



the Markovian full characterization of legitimate 

memory kernels is still missing. Moreover, the concept of 
non Markovianity in this framework is itself not uniquely 



defined 


(jBreuer, Laine and Piilo 


2009 


Rivas, Huelga and 


Plenio 


2010 





The paper is organized as follows: The notation and 
the basics properties of quantum channels are presented 
in Sec. |ll] Section |III| formalizes the notion of memory 
quantum channels and discusses the physical mechanisms 
that are responsible for the arising of correlations in 
quantum communication. The various classes of mem- 
ory channels are reviewed in Sec.|TV] Section [V| discusses 
quantum error correction codes and achievable informa- 
tion transmission rates. Coding theorems and capacities 



are presented in Sec. VI Section VII concerns models 
of memory channels that can be solved for their com- 
munication capacities. Then, quantum channels arising 



from dynamics with memory are discussed in Sec. VIII 



Finally, conclusion and outlook are given in Sec. 



II. QUANTUM CHANNELS: BASIC PROPERTIES 

In a typical communication scenario two parties (Alice 
the sender of the message and Bob the receiver) try to ex- 
change (classical or quantum) information by encoding it 
into (possibly infinitely long) sequences of signals which 



propagate through the medium that separate them. A 
train of transmitted signals defines a sequence of inde- 
pendent uses of the communication line {channel uses), 
and their input-output evolution from Alice to Bob is de- 
termined by the noise which tampers with the transmis- 
sion process. In quantum information the channel uses 
are represented by the degree of freedom (e.g., polariza- 
tion, spins) of a collection Q := {^1,92, • • ■ } of identical 
information carrying objects (e.g., optical pulses, flying 
atoms or ions) which are locally produced by Alice and 
organized in a time-ordered sequence. In this setting the 
noise can then be described by assigning a proper map- 
ping which acts on the (global) input states of the in- 
formation carriers Q to produce the associated (global) 
output states received by Bob, i.e. 



PQ e eiUg) ^ PQ. - $(pq) e 6{nQ, 



(1) 



[here &{T-Lq) and G{T-Lq') stand for the sets of the den- 
sity matrices associated with the Hilbert spaces Hq , Hq' 
which describe the carries as seen by Alice and Bob re- 
spectively - the two need not to coincide in general]. 
Sees. |III| and |IV| will show in which way memory effects 
can arise in these processes focusing on the features of 
the transformation $ that are responsible for introducing 
correlations among the various components of Q. Before 
doing so, however, it is useful to remind that quantum 
mechanics poses some fundamental structural constraints 
on any transformations of the form ([T]), which must apply 
independently from the underlying physical mechanisms 
that is governing the process and independently from the 
composite nature of input system Q. In particular, as- 
suming that prior to the communication, no correlations 
exist between Q and the channel medium, the mapping 
$ must be 

i) linear [i.e. it must transform mixtures of in- 
put density matrices into the corresponding mix- 
ture of the associated outputs, i.e. '^iPiPqii) '-^ 
J2iPi^iPQi''-))j with Pi the probability associated 
with the input state pqii)- This in particular al- 
lows one to expand $ as a linear super-operator on 
the full algebra B{'Hq) of the bounded operators of 

ii) completely positive [i.e. when acting on Q it must 
preserve the positivity of any initial joint density 
matrix of Q with an arbitrary external ancilla 
A that is NOT propagating along the channel]; 

iii) trace-preserving, [i.e. it must preserve the normal- 
ization of all input states]; 

a violation of any of these three conditions implying the 
impossibility of maintaining the statistical interpretation 



of the theory (Breuer and Petruccione 2002 Holevo 



2011 Petz 2008). In the jargon of quantum informa- 
tion theory, any transformations ([I]) that fulfills the three 
conditions given above is said CPTP (completely posi- 
tive and trace preserving) mapping or simply quantum 
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channel. While referring the reader to (Bengtsson and' 


Zycz 


kowskil 2006 


Breuer and Petruccione 


2002;, Holevo, 


2011 


Holevo anc 


Giovannettil 2012| Key 


l| 20021 Petzl 


2008 


) for an exhaustive characterization of their prop- 



erties, one reminds that CFTP maps form a closed set 
under convex combination and a semigroup under super- 
operator composition, [of course when proper correspon- 
dences between the associated input and output spaces 
applies] . It is also well known that any quantum channel 
typically admits unitary dilations that allow one to rep- 
resent it in terms of a unitary coupling with an external 
(possibly fictitious) environment E which not necessarily 
corresponds to the channel medium. For instance, taking 
for simplicity Hq — Hq', one can write 



^{pq) ^ Tte\Uqe{pq ^ \uj)eH)uI 



QE 



(2) 



where \uj)e is a fix state of E, Uqe is the unitary trans- 
formation which couples the latter to the input system 
Q, and Tr^ denotes the partial trace over environment. 
Equation ([2| can also be put in correspondence with the 
operator sum (or Kraus) representation of $, 



d 

$(pq) - V K, PQ k] , 



(3) 



by identifying cIe with the dimension of E and the Kraus 
operator Kj with the linear operator E{i\UQE\^) e of 
Hq, where b; j = 1, . . . , (Ie} is an orthonormal basis 
of E. There exists finally a fundamental relation, known 
as the Choi-Jamiolkowski (CJ) isomorphism fChoi 1975 



Jamiolkowski [1972), which permits to faithfully repre- 



sent any CFTP $ as a density matrix of a composite 
system QA with A being an auxiliary system isomorphic 
to Q. The explicitly connection is obtained by applying 
the map $ to half of a maximally entangled state \I3)qa 
of QA to create the so called CJ state of the channel. 



- ($®id)(|/3)Q^(/3|), 



PQA 



(4) 



[here id stands for the identity map on A\. 

Equation ([I]) implicitly assumes the Schrodinger pic- 
ture in which the states of the system are evolved while 
the observable are kept fixed. In the Heisenberg picture, 
in which instead the states are fixed and the observables 
evolve in time, the channel is represented by the dual 
map O I— $*(0) defined on the space of the receiver 
observable O and identified through the identity 



Tr[<i>(pQ)0] = Tr[pQ$*(0)] 



(5) 



As its Schrodinger counterpart, the transformation $* is 
linear and completely positive, but in general it is not 
trace preserving. On the other hand it is always unital, 
i.e., it maps the identity operator into itself. Operator 
sum representations for $* can be easily constructed from 
those of $ ([S]), yielding. 



$*(0) = gifjOif, 



(6) 



A. Degradable, PPT, and Entanglement Breaking channels 

Associated with the unitary dilation ([2| is the notion of 
the complementary channel of The latter is the CFTP 
mapping $ : &{'Hq) — )■ 6 (He) which moves the initial 
states of the system Q into the states of the environment 
E via the transformation 



PQ ^ $(pq) Trg 



Uqe{pq 8) 



QE 



(7) 



where now Trg denotes the partial trace over the sys- 
tem Hilbert space. This definition allows us to introduce 
another property of quantum channels, which is called 
degradability (Devetak and Shor 2005). Loosely speak- 



ing, a map is degradable when one can recover the final 
environment state ^(pq) just applying a third CFTP 
map to the output system state. More rigorously, a 
degradable map is such that there does exist T CFTP 
satisfying the relation: 



^pq] = iTo<i>){pQ) -.^rmpQ)) , 



(8) 



for any input state pQ [here "o" stands for the super- 
operator composition] . Similarly, a channel is called anti- 
degradable when the opposite relation holds, i.e.. 



*[Pq] = (To^Kpq) :=r(l>(pQ)) . 



(9) 



Another important class of quantum channels is repre- 
sented by the positive partial transpose (PPT) or binding 
maps. These are CFTP transformations <& : &{Hq) — >■ 
&{'Hq) which, when operating on joint states pqa pro- 
duce outputs having positive partial transpose, or equiv- 
alently, which possess a positive partial transpose CJ 
state ( [Horodecki et aL[ |1996[ |Peres[ |1996| |Rains[ |200T |. 
A proper subset of PPT channels is represented by the 
so-called entanglement breaking channels, characterized 



by a separable CJ state ( Horodecki, Shor and Ruskai 
20031 IRuskail [20031). 



B. Causal, semi-causal, and localizable channels 

A couple of interesting notions that have been devel- 
oped in Refs. ([Beckman et aL[[2b01 Eggehng, Schlinge- 



mann and Werner) 2001 Piani et aL[ |2006) for the special 
case in which Q can be treated as a finite bipartite sys- 
tem, e.g. by splitting the set of carriers {qi,q2,- ■ ■ } into 
two independent subgroups Qi, Q2 so that Q = QiQ2- 
In this context the transformation $ of Eq. ([l]) i s said 
to be Qi Q2 semicausal (Beckman et al. 2001) if for 



any local quantum channel 4* applied to Qi before the 
action of the map <&, there is no detectable effect in the 
subsystem Q2, i.e. 

Trg, [$ (pq)] = Trg, {$ [(* ® id) (pg)]} , (10) 

where pg is a generic [possible entangled] input state 
of the two carriers and where Trg^ [• • • ] describes the 
partial trace with respect to Qi. In other words, for 
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Qi ~^ Q2 semicausal maps cross-talking from Qi to Q2 
is prevented. Similarly, one may introduce the notion of 
Q2 ^ Qi semicausal map. When both properties are 
satisfied, the map is called causal or non-signaling. Spe- 
cial examples of non-signaling channels are the tensor 
product channels which can be expressed as $ = $1 (g) $2 
with $12 being CPTP maps operating locally on Qi and 
Q2 respectively. 

Still remaining in the bipartite scenario discussed 
above ^ is said to be localizable if it can be implemented 
by means of a bipartite state shared between two par- 
ties that operates locally on the sybsystems Qi and Q2 



without communicating (Beckman et al. 2001). That is. 



localizable channels do not need communication between 
Qi and Q2 to be simulated. Formally, this is the case 
when $ satisfies the following property 

*(PQiQ2) =TrAiA, [(*®f^)(pQiQ. , (11) 

where (JJA1A2 is a shared bipartite state, and ^ and f2 
are now quantum channels acting locally on subsystems 
QiAi and Q2A2, respectively. Otherwise, $ is called 
Qi — ^ Q2 semilocalizable if one-way communication from 
Qi to Q2 is required. Accordingly in this case Eq. (11) 
is replaced by 



^(PQiQ2) = Tr^ [(* ° ^){pQiQ2 ® ^a)] 



(12) 



with uja being the state of an ancillary system A that 
acts as the mediator between Qi and Q2, while ^ and 
are quantum channels acting on the systems Q2A and 
QiA respectively. Analogously, one defines a Q2 — > Qi 
semilocalizable map. 



By comparison of (11) and (12) it follows that all 
Qi ^ Q2 semilocalizable maps are Q2 ^ Qi semicausal, 
which in turn implies that all localizable maps are causal. 
Moreover, it can be proven that semicausality implies 
semilocalizability, hence semicausal and semilocalizable 
maps coincide, although causal and localizable maps do 



not (Beckman et al. 2001 Eggeling, Schlingemann and 


Werner 2001 Piani et al. 


2006p. 



each element of the sequence Q. Under this assumption 
the multi-use map associated with the communication 
line is expressed as a tensor product of a CPTP map 
$ : QCHq) ^-> QCHq) that acts on the states of a single 

carrier q. Therefore indicating by H^q^ := Hf"' — Hq^ (8) 
■■■(E) Hq^ the Hilbert space of the first n carriers of the 

system, its input density matrices pg*^ G ©('Hq''') will be 
mapped into 



$(«)(>)) ^ 



(13) 



with := <i>(g) • • • (g) $. Equivalently, we can say that 
the Kraus operators of the memoryless map can be 
expressed as a tensor products Ki^ (g) • • • (g) Ki^ formed 
by independent and identically distributed sequences ex- 
tracted from the Kraus set {Ki}i associated with the sin- 
gle carrier channel $. For an explicit example consider 
for instance the model discussed in Ref. (Giovannetti 



2005). Here the carriers Q are assumed to propagate 
from Alice to Bob, one by one and at constant speed, 
while interacting with an external environmental system 
via a constant coupling described by the unitary operator 
Uqe £ B{%q (g He). The environment e is also assumed 
to undergo a dissipative process which, on a time-scale 
T tends to reset it into a stable configuration uje- The 
memoryless regime is achieved in the limit in which the 
rate v at which the carriers propagate from Alice to Bob, 
is much lower than the inverse of the relaxation time 
T, i.e., v <C 1/t. In this case in fact each carrier cou- 
ical environmental states. Defining then 
• • g) uje^ , this allows us to write 

^^9„,e„ g) ■ • • g) f/gi,ei 

{p^^^^u;T)Ul,e.®---^Ul^J, (14) 



UJ 



On 



pj^ ^ Tr^ 



which reduces to Eq. (T3| when identifying Ti e[Uqe{^ ■ ■ (g 
We)C^ge] with unitary dilation of the single-use channel $. 

B. Non-anticipatory memory quantum channels 



III. FROM MEMORYLESS TO MEMORY QUANTUM 
CHANNELS 

Having in mind the multi-uses communication scheme 
detailed at the beginning of Sec. [H] in which a time- 
ordered sequence of carriers Q := {gi, 92, • • • } propagates 
from Alice to Bob along a noisy channel, in this section 
one starts reviewing the mechanisms that are responsi- 
ble for the introduction of correlations [memory effects] 
among the various elements of Q. 



A. Memoryless quantum channels 

Memoryless quantum channels describe those scenarios 
in which the noise acts identically and independently on 



Whenever the tensorial decomposition of Eq. ( 13 ) 
doesn't apply, one can speak of memory channels or 
correlated noise channels. Among the plethora of pos- 
sibilities, the following will focus only on those config- 
urations which have physical relevance and which have 
attracted some interest in the recent literature. In par- 
ticular, one shall treat those models in which the noise 
respects the time-ordering of the carriers Q so that at 
a given channel use, the output cannot be influenced by 
successive inputs. Under this condition there must exist 
a family of CPTP maps J" = 1, 2, • • • } with 

: Gin'-Q^) &{'H'-q^) which allows us to express 
the output states of the first n carriers in terms of the 
density matrices of their associated inputs, i.e., 



(15) 
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Clearly the property ( 15 ) requires that the family T must 



m^, and write Uq-u as a transfor- 



fulfill the minimal consistency requirement that for all 
m < n the element should be obtained as a restric- 
tion of 'f'^"'^ over the degree of freedom of the first m car- 
riers. That is, given p^"^ e ©(H^"^) and pj"' e &{'H'q^) 
its reduced density matrix associated with the first m 
element of Q, one must have 



terms \Q)rm |0)m„, ^^^^ ..^i.^^ ^q^- 

mation which couples the j-th carrier only with the first 
j elements of M, i.e., 



(18) 



$(™)(>^)) 



(16) 



with Im' being the identity operator on the m' compo- 
nents of the environment, see Fig.[I]b). 

An alternative, but fully equivalent, representation 
for non-anticipatory channels is obtained by adding to 



where Tr'™-* stands for the partial trace over all the car- 
riers but the first m. Borrowing from the classical the- 
ory of communication ( |Gallager[ |i968| one can dub these 
communication lines non- anticipatory quantum channels. 
In Ref. [Kretschmann and Werner[[2005| non-anticipatory 
channels are also called causal, it is worth stressing how- 
ever that this definition of causality has nothing to do 
with the one introduced in Sec. III.BI 



1. Structure theorems for non-anticipatory quantum channels 



Eq. (|17| a collection of local environments which indi- 

(19) 



vidually couples with the carriers, i.e. 



$(")(p(j')) 



-Trji/B 
[Pq 



u, 



qnMe„ 



■u, 



giMei 



■u. 



q„AIe„ 



where for all j — 1,2, 



U, 



qjMe. 



is now the unitary 



transformation which describes the coupling of the j-th 
carrier with its own local environment Cj and with the 
memory system M, where w^" := (8) • • ■ We„ as 
in the memoryless case, and Qm is some fixed state of 
Af, see Fig, [l] c). Equation (p]) was first introduced 



In the language introduced in Sec. |II.B| non- 
anticipatory channels can be classified as semicausal 
with respect to the natural ordering of the channels 
uses. A direct generalization of t he semicausal/ semilocal 



connection discussed in Refs. (Beckman et al. 



Eggcling, Schlingemann and Werner 2001 Kretschmann 



2001 



in Ref. |Bowen and Mancini[ |2004| as a model for rep- 
resenting correlated channels: from Eq. (17) it follows 
that it provides a general unitary dilation for every non- 
anticipatory quantum maps. It can also be expressed in 
terms of a n-fold concatenation of a sequence of CPTP 
maps acting on a single carrier and the memory system 



and Wernerl 120051 IPiani et a/.l 120061) can then be used to M (|Bowen and Mancini| |2004| [Kretschmann and Werner] 



show that these channels admits a semilocal representa- 
tion where each carrier couples sequentially with a com- 
mon memory system M, whose back-action on the mes- 
sage state simulates the memory effects of the transmis- 
sion. Accordingly, all the non-anticipatory CPTP maps 
can be expressed as 



2005). Such concatenation is shown pictorially in Fig. 11 



c) and results in the following identity 



I 



$(n)(>))^TrM 



^^QlApi'^^^M) , (20) 



with 



<i>(")(p(")) = Tr 



M 



Uq^M ■ ■ ■ Uq^M 
in) 



(») 



o $ 



(n-l) 
QM 



^q^M O ... O ^q^M , (21) 



{p^Q^^nM)ulM---ulM\ : (17) 

where for all j = 1, 2, • • • , n, Uq^M is a unitary transfor- 
mation which describes the coupling of the j-th carrier 
with the memory system Af , and where J7m is some fixed 
state of Af, see Fig. [I] a) . An explicit proof of E q. (17) 
was first given in Ref. [Kretschmann and Werner| |2005| in 
the context of quasilocal algebras, under the assumption 
of translational invariance of the noise [more on this will 
be provided in Sec 



where for j = 1, . . . , n, $ 



93 



M 



G{Hq^(E}H 



LM 



eiUq 



T-Lm) is a CPTP map that operates on the j-th carrier 
and on the memory ancilla Af and which is defined by 
the unitary dilation 



Tr, 



Uq^Me,{- 



iMe 



(22) 



III.C 



An alternative proof can also 

be found in Appendi 

In Eq. (17) Af is in general a large system whose di- 



mension d]\i is an explicit function of n (in any case it 
can always be chosen to be less than or equal to d^" with 
d = dimHq being the dimension of a single carrier). As 
a matter of fact, as explained in Appendix [C] one can 
take Af to be a composite system of components mi, 
, • • • , whose dimensions can always be chosen to 
be not larger than d^ . In this configuration then one can 
assume ^Im to be a pure tensor product state of local 



Cases of special interest ( Kretschmann and Werner[ 
2005 ) are those in which, for all j, the ^q-M describes the 
same mapping <I> = on &{%q®'HM) which, accord- 
ing to Eq. (21 ) becomes the generator of the n-fold con- 
catenation. Memoryless channels can then be included in 
this class as limiting case in which the generator $ can 
be expressed as a tensor product channel that acts inde- 
pendently on the carrier q and on the memory system Af. 
In terms of the unitary dilation ( 19 ) this is equivalent to 



assume that the unitaries Uq-Mcj in Eq. (14) factorize in 



a tensor product t/(8) Vm, where Vm is a unitary operator 
on the memory system and J7 is a unitary acting only on 
the degree of freedom of the j-th carrier and on its local 
environment e,-. 
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FIG. 1 Unitary dilations for a non-anticipatory quantum 
memory channels, a) graphical sketch of the representations 
of Eq. (17 1: here the noise correlations among the n chan- 



nel uses can be described via a series of concatenated unitary 
interactions with a common reservoir M whose dimension in 
general depends (exponentially) upon n (n = 3 in the ex- 
ample). Notice that while the carrier gi might influence the 
outcome of 5293 via their common interaction with M, 5293 
cannot influence the output of the first carrier; b) the envi- 
ronment M can be also represented as a collection of smaller 
systems Mi , M2 , • • • initially prepared into a separable state 
while, as shown in Eq. (18 1, the unitary transformation op- 
erating on the j-th channel use couples it with the first j 
subsystems only; c) unitary dilation ( |19[ ) where apart form 
M a series of local environment ei, 62, ■ • ■ , are also present. 
In all the diagrams the unitary operators (represented by the 
white boxes) are applied sequentially on the input states of the 
global system (i.e., the carriers and the environment) starting 
for the one on the top of the figure. The carriers and the en- 
vironmental states evolve, respectively, from-left-to-right and 
from-top-to-bottom while interacting meeting at a white box. 
The trash-bin symbol stands for the partial trace operation 
on the environment. 



2. Finite-memory channels 



The word finite-memory (Bowen and Mancini 2004) 



has been used to indicate those non-anticipatory channels 
which admit a representation of the form ( 19 ) with M be- 
ing finite dimensional. The dimension of the memory is 
determined by the number of Kraus operators in the sin- 
gle channel expansion and the correlation length of the 
channel, which may be defined as the maximum number 
of channel uses for which the noise is not conditionally in- 
dependent. Any channel with a finite correlation length 
may be generated by a channel with a finite memory, ac- 
cording to this model. Within the representation (17) 



examples of finite-memory channels are obtained by as- 



suming that the unitary transformations ( 18 1 couple the 
carriers with no more than a fixed number k of environ- 
mental subsystems, the parameter k playing the role of 
the correlation length of the channel. More precisely for 
all j > k one has. 



(23) 



Ql 



92 



^3 



O Q 



94 



for a finite-memory non- 
channel with correlations 



FIG. 2 Unitary dilation (|17l 
anticipatory quantum memory 
length k = 2. In the depicted example the total channel 
uses is n = 4 and each carrier is supposed to interact with 
only two components of the environment. Consequently the 
carrier gi can influence the output carrier ga only via g2 . For a 
comparison see the scheme of Fig.[l]b) where instead the flrst 
carrier can directly influence ga via their common interaction 
with mi). 



with (see Fig. [2] for a graphical representation of the case 
with k = 2). Extreme case of finite- memory channels is 
provided by the memoryless one where k — I and each 
carrier interacts with a devoted component of the multi- 
partite environment M (specifically, for each j, the car- 
rier Qj interacts with rnj only). 



3. Symbol independent vs intersymbol interference maps 

A special subset of the non-anticipatory channels is 
formed by the symbol independent (SI) maps. They are 
communication lines where previous input states do not 
affect the action of the channel on the current input state. 
In other words the symbol independent maps are non- 
anticipatory (or semicausal) with respect to all possible 
ordering of the carriers. Therefore given a generic sub- 
sets of the carrier set Q, its output state is uniquely de- 
termined by the corresponding input state via a proper 
CPTP mapping. Following the notation introduced in 



Refs. ([B cckma n et aL||2001 



Werneri ,2001; iPiani et al. 



Eggeling, Schlingemann and 



2006), they can be said non- 



signaling (or causal) channels, meaning that the output 
states of any subset of the carriers cannot be influenced 
by the input state of the remaining carriers. 
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This is opposed to the channels with intersymbol inter- 
ference (ISI), where the input states of previous carriers 
affect the action of the channel on the current input. An 
extreme example of such channels is the quantum shift 
channel, where each input state is replaced by the previ- 
ous input state. 



C. Quasi-local algebras approach 

Till now one has followed a constructive approach 
to quantum channels with memory, where CPTP maps 
which process long messages were always thought of 
as concatenations of smaller units which, starting from 
an official "first carrier" element, process one quan- 
tum signal each. An alternative view is proposed 
in |Kretschmann and Werner} |2005| by treating the chan- 
nels as mappings applied on an infinitely long message 
strings. This approach requires some advanced mathe- 
matical tools, like quasi-local algebras, that are sketched 
in Appendix [D] 



1. The Structure of Causal Channels 

To set the stage, suppose to have a quantum chan- 
nel which transforms input states of an infinitely ex- 
tended quantum lattice system [representing the infi- 
nite mess age string] into output states on t he same sys- 
tem. In (Kretschmann and Werner 2005) this map is 



formally assigned by working in the Heisenberg picture 
via the introduction of a completely positive and uni- 
tal map — > operating on the quasi-local alge- 
bra that defines the observable quantities on the lat- 
tice. In this context one says that the channel is trans- 
lational invariant if commutes with the shift oper- 
ator on the lattice. Requiring then that future inputs 
should not affect past measurements [non-anticipatory 



constraint (15)], Ref. (Kretschmann and Werner 2005) 



introduces the definition of causal channel as a com- 
pletely positive and unital translational invariant map 
<I>* that fulfills the constraint 

^* ( l^i-oo^z] ^ j[l + 3,oo)\ ( l,{-oo,z]\ ^ ill + z,oo) 



(24) 

for all z e Z and e Example of 

causal maps are provided by concatenated memory chan- 
nels { Kretschmann and Werner 2005 ) which can be eas- 



ily constructed by adapting the concatenation scheme of 
Eqs. (20)-(21| to the quantum lattice formalism. Inter- 



2. Ergodic channels with decaying input memory 

The notion of channels with memory having a nonan- 
ticipatory and decaying character can be extended to 
quantum channels having ergodic properties. Actually 
such properties have been established till now(Bjelakovic 



and Boche 2008 1 only for classical-quantum channels, i.e. 



CPTP maps that converts classical inputs into states of 
a quantum system ( [Holevo and Giovannetti[ 2012 1. 

To give an idea of how things work for classical- 
quantum channels (cq-channels), let A be a finite set and 
let T-l he a d-dimensional Hilbcrt space. By one de- 
notes the set of doubly infinite sequences with compo- 
nents from A and is the quasi-local C* -algebra cor- 
responding to the set of bounded operators B on H (see 
Appendix |d]) . Then, a cq-channel is usually meant as a 
triple {A, Ti, W), where W maps elements of A into den- 
sity operators on "H. By the quasi-local algebra formalism 
it can be regarded as the map W : A^ x B^ ^ C such 
that for each 6 G B^, W{x,b) is measurable and for each 
X £ A^, W{x, b) is a state. 

A cq-channel W is said to be stationary if W{TinX, b) — 
W{x,Tontb) holds true for all x & A^ and all b e B^. 
Here Ti„ and Tout denote the shift on the input and out- 
put algebra respectively. A cq-channel W is called er- 
godic, if it is extremal in the convex set of stationary 
cq-channels. As such it cannot be written as a statis- 
tical mixture of cq-channels W = J^iPi^i- Following 
the definition of causality given in Section |III.C.1[ a cq- 
channel W is causal if for each n £ Z, b e and 
all x,x A^ [xi = Xi'jVi <n), it is 



Wix,b) = W{x,b) . 



(25) 



Furthermore, W is called input memoryless if for each 
neZ,be and ah x,£eA^ (x, = i^; Vi > n) it is 



W{x,b) = W{S:,b) . 



(26) 



In the presence of memory it s expected that the effect 
of the far past inputs do not affect much present and 
future outputs. According to this idea, a channel W is 
said to have decaying input memory (and anticipation) if 
for each e > there exist non-negative integers m(e), a(e) 
such that 

\W{x,b)-W{x',b)\<e, (27) 

for all b e ^l^'^+fe], n,e Z,k e N, whenever Xi = x'^ 
for n — TO < i < n -\- k -\- a and m > m(e), a > a(e). 
This provides a 'continuity' property of the channel which 
plays a crucial role in establishing coding theorems, an 
idea that will also appear in Section IV. D and goes back 
to the classic paper by McMillan (iMcMillanl |l953l). 



estingly enough Ref. (Kretschmann and Werner 2005) IV. TAXONOMY OF QUANTUM MEMORY CHANNELS 



proves a structure theorem similar to the one discussed 
in Sec. III.B.l which shows that any causal map (24) can 



always be represented as concatenated memory channel 
with a uniform generator. 



This section focuses on some classes of non- 
anticipatory memory quantum channels which have at- 
tracted special interest in the literature. 
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A. Localizable, fully non-anticipatory memory quantum 
channels 



many-body environment 



A subset of 



quantum channels 



which have been extensively analyzed (Giovannetti and 


Mancini 


2005 


Plenio and Virmani 


20071 2008|) is graph- 



ping ( 15 ) can be expressed in terms of (not necessarily 
identical) local unitary couplings with a correlated many- 
body environmental system E := {ei, 62, • • ■ }. Memory 
effects in the communication appear because, differently 



from the memory less case (14), the many-body environ- 
ment is not initialized in a state oj^-* which does not 
factorize, i.e., 



$(")(p(") 



Tr, 



(n) „ (71] 



(28) 



An alternative representation of these maps has been 
provided in Ref. (Caruso, Giovannetti and Palma 2010) 



by generalizing a model presented in Ref. (Ban, Sasaki 



and Takeoka 2002 Bowen and Bose 2004) for memo 



ryless channels. According to this approach the channel 
noise is ef fectively described as a quantum teleportation 
proto col ([Bennett et al. 1993 Braunstein and Kimble 



1998 Vaidman 1994) that went wrong because the com- 
municating parties used non optimal resources (e.g., the 
state they shared was not maximally entangled). In the 
case of ( 28 ) each of the carriers gets teleported indepen- 



dently using the same procedure, the correlations aris- 
ing from the fact that the communicating parties use as 
shared resource a correlated many-body quantum state. 

One can easily verify that the maps described by 
Eq. (28 1 are symbol independent. More precisely in 



the language of Refs. (Beckman et al. 




'2001 Eggcling, 


Schlingemann and Werner, j2001 


Piani et al. 


.2006) they 



are localizable, non-signaling channels, since by construe 
tion admit a representation in terms of local coupling 
with a correlated ancilla (it is worth reminding that not 
all the non-signaling maps can be represented in this 
form, see e.g. (Beckman et al. 2001 Piani et al. 2006 



Popescu and Rohrlich 20011)7 



B. Perfect memory channels 

The memoryless channels are characterized by pos- 
sessing unitary dilations in which the environment has 
a dimension which is exponentially growing in n (i.e., 

log[dimH^''] = nlog[de]) or, equivalently, by possessing 
a (minimal) operator sum representations whose Kraus 
sets contains a number of elements which is exponen- 
tially growing in n. This same property typically holds 
also for memory channels with the important exception of 



the perfect memory channels ( 


Giovannetti, Burgarth and 


Mancini 


2009 


Kretschmann and Werner 


2005 


). The 

Dtained 


simplest example of such commimication lines is 



by assuming that the memory system M in Eq. (17) does 
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FIG. 3 Model for a localizable, fully non-anticipatory quan- 
tum memory channel. Here the correlations are introduced 
by allowing the state of the environment (gray element) to be 
initially entangled. As in the previous figures white boxes rep- 
resents unitary couplings while the trash-bin indicates partial 
trace over the corresponding degree of freedom. These maps 
are SI and hence non-anticipatory (therefore they also admit 
unitary dilations of the form described in Fig. [TJ. 



not scale with n and it is finite dimensional. Under this 
hypothesis the maps 4>*^"-' explicitly admit a unitary di- 
lation with an environment (the system M) of constant 
size which is hence unable to capture all the information 
that is sent through the channel (the latter being mea- 
sured by the size of the carriers "H^^ which is exponential 
in n). As a consequence, in the asymptotic limit of long 
carrier sequences, no information is expected to be lost 
to the environment, yielding optimal quantum commu- 
nication efficiency. A typical example is provided by the 
shift channel (see Sec. 



III.B.3) which can be described 



as in Eq. (17) by assuming M to have the same dimen- 

as the SWAP 



M 



sion of a sing 

operator. It is worth noticing that in Ref. (IBowen, De- 



e carrier and by taking Uq 



vetak and Mancini 2005) it was also conjecture that the 
types of memory channels that, analogously to the shift 
channel, display only intersymbol interference may be 
represented as perfect memory channels. 

More generally the class of perfect memory channels 
can be extended to include all the CPTP maps (15) that 
admit unitary dilations (17) in which the dimension d-M 
of the environmental system M is sub-exponential in n, 
i.e., 



lim -log[c?M] 

n— J-oo n 



(29) 
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As discussed explicitly in Sec. VLB. 2 also in this case the p' 



channel capacity of the associated communication line 



is optimal ( 


Giovannetti, Burgarth and Mancini 


Kretschmann and Werner 2005 ) . 



C. Markovian channels 

An important class of non-anticipatory quantum chan- 
nels is given by the channels with Markovian correlated 
noise. They describe noise models in which the carri- 
ers are transformed via the applications of strings of lo- 
cal CPTP maps whose elements are randomly generated 
by a classical Markov process. For instance a quantum 
Markov process is given by 

x$('.)o$(^-i)o...o$(;^)(pf)),(30) 

where is a set of CPTP maps operating on the 

j-th carrier, where is an initial probability distribu- 
tion, and where finally for j >2, the p'f^l are conditional 
probabilities. 

The mapping ( 30 1 is symbol independent since mod- 
ifying the input state of previous channels uses do not 
have any effects on the output states of the carriers that 
follows. A unitary dilation of the form ( 19 1 can be ob- 
tained by identifying the initial state m of the memory 

M with the pure vector J2i \fp^\'^)M, and by taking the 
unitary Uq-Me^ in such a way that for all vectors IV')?^ of 
the j-ili carriers one has 

[/,,MeJV',*',0) = 4f , (31) 

for J = 1, while 

i7,^.Me,#,*',0) = (32) 

i 

for j > 2 (in the above expressions {A^q^{l)}i is 
a Kraus set for i, i^i/^^) stands for the state 

\ij;)q.\i)M\4'%',de,, while {\i)M}k and {\(j)i,i)ej}i,i are or- 
thonormal basis for the memory systems M and Cj re- 
spectively). 

Most of the analysis conducted so far focused on the 
special case of homogeneous Markov processes in which 
both the and the <I>g/ do not depend upon the car- 



rier label j (e.g., p^^, := Pi\i'). Under these conditions 
one also says that the quantum Markov process is reg- 
ular if the corresponding classical Markov process p^^i 
is regular, i.e., if some power of the transition matrix 
r of elements Pi\i' has only strictly positive elements. 
In this case, for j ^ oo the statistical distribution of 
the local noise converges to a stationary distribution 



lim 



«^oop1"^ with 

(«) STr-nn-i-, (1) 



Pi 



(33) 



being the probability of getting ^q^ 



Also the initial probability pl 



(1) 



on the j-th carrier, 
is said to be stationary 



if it satisfies the eigenvector equation J2i> ^iyPi/ ^ = Pi^^ 



(when this happens pf^ = pl'^' and the local statistical 
distribution of ^q^^ is identical for all the carriers). 

The first example of a regular Markov process has been 
analyzed by |Macchiavello and Palma,,2002 , Here the car- 
riers are assumed to be qubits and the CPTP transforma- 



(1) 



tions $q entering in Eq. ( 30 ) perform unitary rotations 
. . ^Yiere af^ = I is the identity 
operator while for i = 1,2,3, cTg''' is the i-th Pauli spin 
matrix. The conditional probability Pi\ii which describes 
the associated classical Markov process was finally writ- 
ten as Pi\i/ = (1 — m)Pj-^^ + tJ-Sii' where fi G [0, 1] is a cor- 
relation parameter (notice that for /z = the model de- 
scribes a memoryless channel while for = 1 it describes 
a long- memory channel - see Sec IV.E). This model of 



Markovian correlated Pauli channel show a remarkable 
feature when it is used for the transmission of classical 
information (see Sec. VI). That is, when two successive 
uses of the channel are considered, classical information 
is optimally encoded in either separable states or max- 
imally entangled states, depending whether the correla- 
tion parameter fi is below or above a certain threshold 
value. This feature was first conjectured in ( Macchiavello 



and Palma 2002), then proven for certain instances of 



the model in (Macchiavello, Palma and Virmani |2004 1, 



and finally proven for general Markovian correlated Pauli 



channel in (Daems 2007). 



An experimental demonstration of the optimality of 
entangled qubit pairs for encoding classical information 
through a correlated Pauli channel was provided by |Ba- 



naszek et al. 2004, for mechanically induced correlated 
birefringence fluctuations, which in turn induce corre- 



lated depolarization ( Ball, Dragan and Banaszek 2004 1 . 



A generalized model of d-dimensional Markovian cor- 
related Pauli channel was considered by |Shadman et al. 



|2011| for the problem of sending classical information us- 
ing a dense-coding protocol. An alternative model of two- 
qubit correlated channel was characterized by |Caruso ei\ 
2008b| in terms of the minimum output entropy. 



al 



Going beyond the case of two uses of a qubit chan- 
nel, Markovian correlated depolarization over an arbi- 



trary number of channel uses was studied in ( Demkowicz- 
Dobrzanski, Kolenderski and Banaszek[ 2007 Karim- 



ipour and Memarzadeh 2006b), and the case of Marko- 



vian correlated noise in higher dimensional quantum sys- 



was considered in ( 


Karimipour and Memarzadeh 


Karpov, Daems and Cerf 


2006|b 


) . Generally 



speaking, the optimality of entangled state encoding can 
be interpreted as a consequence of a decoherence-free 
subspace (see Sec. |V| associated to the correlated noise 
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model: this has been considered for the Hilbert space de- 



fined by multiple uses of a qubit channel in ( 


Demkowicz- 


Dobrzanski, Kolenderski and Banaszek 


2007 


) and for the 



multiphoton Hilbert space associated to the polarization 
of light ( jBall and Banaszek, ,2005^ . In a different con- 
text, the same phenomenon has been discussed for the 
problem of quantum communication with polarized light 
without a shared reference frame (Bartlett, Rudolph and 



Spekkens 2003 1 



Finally, models of Markovian correlated noise in the 
framework of quantum systems with continuous vari- 



able (seeSec. VII. B I were first discussed in ( Cerf et al.^ 
20051 |2006l ) for the case of two uses of the channel, and 



then extended to the arbitrary number of uses in (Lupo, 
Memarzadeh and M ancinT, |2009| [Schafer, Karpov and 
CerfH2009[) (seeSec. [vIlb|). 



D. Fixed-point, indecomposable, and forgetful channels 

Memory channel representations where the action on 
the memory state is independent of the input state are 
termed fixed-point memory channels, as the memory 
state will hav e a fixed point under the action of t he repre- 
sentation ( ,Bowen, Devetak and Mancini 2005). Specif- 
ically remembering the definitions of Eqs. (20) and (21) 
this implies, 



nM{p,n) Trg ® ^m) 



M 



(34) 



for all n and for all Pq ■ Fixed-point channels can be 
easily shown to be symbol independent while the op- 
posite is not necessarily true. Indeed from Eq. (21) 
one has that the output state of n-th carrier p^^^ :— 

TrQ(„-i) <i>(")(p[?-') can be expressed as 



Tro(„-i 



Tr 



M 



(35) 



which only depends upon the reduced density matrix 
(in these expression TTQ(n-i)j^,j[- ■ ■] indicates the partial 
trace with respect to M and the first (n — 1) carriers). 

An indecomposable channel is one where, for each chan- 
nel input, the long-term behavior of the channel is inde- 
pendent of the initial memory state ( Bowen, Devetak and 
Mancini 2005[ ) . More precisely a finite-memory quantum 
channel is indecomposable if for any input state p and 
e > there exists an N{e) such that for n > N{e), 



D{QM{p,n),T,M{p,n)) < e , 



(36) 



where flMip,n) and Sm(/3, "-) are the memory states af- 
ter n uses of the channel for the intial memory states 
rijii and Sjv/ respectively, and where D{p, p') is the trace 
distance (see Appendix [A| . 



The main features of indecomposable channels have 
been revisited through the notion of forgetful chan- 
nels ( [Kretschmann and Werner , 2005 1 which it is pre- 
sented here in a slightly more general form. Let (^qM be 
a CPTP map on &{Hq (E) Hm) which generates a non- 
anticipatory channel ^qI.j via the n-fold concatenation, 
Eqs. (20l-(22). We say that ^qM is forgetful iff there ex- 
ists a sequence of quantum channels Sn ■ &{'Hq'SS'HM) 
GCHm) such that 



lim ||<i> 



QM 



in) 



- iAm ® SI 



= 



(37) 



where 



is the cb-norm (see Appendix IaI) while 



*^QA/^"'' Sn are the Heisenberg dual of the maps <&q2/ 
and Sn respectively. An example is given by the channel 



^^q\.j obtained via concatenation (22) using as generator 
the map <^qM := p id -|- (1 -p) SWAP, where p e [0, 1), 
and SWAP denotes the swap channel which exchanges 

(n) 

q and M. The only way for ^qIj not to be forgetful is 
to choose the ideal channel in every step. However, the 
probability for this event vanishes in the limit n — > oo as 
p", implying that Eq. (37) holds. 



The equivalence between the notion of indecompos- 
able channels and that of forgetful channels has been 
shown in ( Kretschmann and Werner 2005 1 . Several cri- 



teria for a quantum memory channel to be forgetful ex- 
ist [for instance, a sufficient condition is to have that 

ll*^QAf ~ i^Af ® S*j\\cb falls below 1 for some finite n]. 

From a physical point of view, one could expect a 
generic quantum memory channel to be forgetful. In- 
deed, it can be proven that the subset of forgetful chan- 
nels is dense and open (according to th e topology induced 
by the norm of complete boundness) ( Kretschmann and 



Werner 2005 ) . In the case of Markovian channels, the 
forgetfulness is determined by the asymptotic properties 
of the underlying Markov chain: in particular, for a dis- 
crete variable memory system, the channel is forgetful if 
and only if the underlying Markov chain converges to a 
unique stationary state ( Datta and Dorlas 2009 ) . On the 



other hand, if the memory system is described by con- 
tinuous variables, one could have situations in which the 
Markov chain has a unique stationary state, yet the con- 
vergence property (37) is not verified. To overcome this 



limitation, a weaker notion of forgetfulness, named weak 
forgetfulness, has been introduced in ( Lupo, Memarzadeh 



and Mancini 2009). This notion applies only to fixed- 
point memory channels with memory systems described 
by classical degrees of freedom. Although restricted to 
this setting, its definition coincides with that of inde- 
composability [Eq. (|36|], and is equivalent to forgetful- 
ness for discrete variable Markov chain. It has also been 
argued in ( |37[ ), but not proven rigorously, that the cod- 
ing theorem for forgetful channels (see Sec. VLB. 4) can 
be extended to weak forgetful ones. Beyond this setting. 



the model of Gaussian mem ory channel in (Lupo, Gio- 
vannetti and Mancini 2010) was proven to be forgetful 
under conditions on the memory initialization, e.g., if the 
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initial state of the memory is a Gaussian state with finite 
first and second moments. Finally, the relation between 
the forgetfulness of the channel and the chaotic quan- 
tum evolution of its environment was studied in ( Barreto 



Lemos and Benenti 2010) for the case of a model of de- 



phasing channel with memory. 



E. Long term memory channels 



Long term quantum memory channels describe those 
communication lines in which the effect of the memory 
does not die with the channel uses. Extreme examples 
are provided by statistical mixture of memoryless chan- 
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I where the carriers are transformed accord- 



he mapping 



(38) 



where is a set of single carrier CPTP channels, and 

Pi is a probability distribution on such set. Equation ( 38 ) 



describes a scenario in which, with some probability pi 
all the carriers of the system are operated by the same 
identical local transformation The index i can be 
interpreted as a "switch" selecting different memoryless 
channels, and the memory channel in (38) can be inter- 



preted as the average channel over different values of the 
"switch" . 



It is worth stressing that the transformation (381 is 
fully non-anticipatory (i.e., the input state of any sub- 
set of carriers cannot influence the output state of the 
remaining ones): as a consequence, fixing an ordering, 
it can always be represented as in Fig. [l] with a proper 
choice of the unitary couplings. Also the memory corre- 
lations of this class of channels can be considered to be 
given by a Markov chain which is aperiodic but not irre- 
ducible ( Norris 1997 ) . This can be easily seen by notic- 

J L (i) 

ing that Eq. ( 30 ) reduces to Eq. ( 38 1 by assuming the 
appearing in the former expression not to depend upon 
the carrier and by setting p'f^^, = Sui for all j — 2, • • • n. 
Hence, once a particular branch, i — 1,...,M, has 
been chosen, the successive inputs are sent through this 
branch: transition between the different branches (which 
correspond to the different states of the Markov chain) is 
not permitted. 



V. QUANTUM CODES 

The channel formalism introduced in the previous sec- 
tions provides the tools for describing the effects of noise 
on a quantum system. Within this formalism the Kraus 
operators have to be intended as errors introduced onto 
the message represented by the system's (input) state. 
Given that one aims at transmitting (classical or quan- 
tum) information through the quantum channel, the 



presence of noise may impose severe limitations to the 
reliable information transmission. The key idea, mutu- 
ated from classical information theory, to realize reliable 
information transmission in the presence of noise is to ex- 
ploit redundancy by means of a suitable quantum error 
correcting code (QECC). The latter must be designed to 
make the received noisy signal decodable with an asymp- 
totically negligible probability of error. 

While classical works on QECC were mainly concerned 
in counteracting noise arising from memoryless quan- 
tum channels, the effects of noise correlations and mem- 
ory characterizing quantum memory channels posed new 
challenges for the optimal design of QECC. 

More specifically, one says that information can be reli- 
ably sent through a quantum channel : ©(C^*^") — > 
©(C^**") if there exist an encoding map 



: 6( 



with n > k, and a decoding map 



©( 



: ©(C^S"') ^ 6( 



such that V|-0) G 971 it is 



(39) 



(40) 



1 



lim F 1 |<^/yv'i,-±-D 

(41) 

Here 9JI is the message space corresponding to 
{|0), |1)}'*'^ in case of classical information transmission, 
and to C^***^ in case of quantum information transmis- 
sion. Furthermore, in Eq. (41 1 F is the fidelity (Jozsa 
1994[ |Uhlma^|1976[ ), defined as in Eq. ([A3|. With the 
above conditions fulfilled, one says that R ~ k/n is a 
reliable communication rate. We will address in the next 
Section the question of what are the maximum reliable 
rates for a given quantum channel. Here, one is inter- 
ested on the image of the encoding map denoted 



C C C which is called quantum error correcting code 
(QECC). The state of k qubits is encoded in C and has 
to be decoded from C with asymptotically unit fidelity. 



A. Standard quantum coding theory 

Let us consider a quantum memoryless channel 
characterized by Kraus operators {^i}, on the Hilbert 
space C^®" of n qubits 



(42) 



where Ai = Ai-^ (E) . . . (E) Ai^ describes independent and 
identically distributed errors on single qubits. A [[n. A:, d]] 
quantum correcting code C is given by a 2''-dimensional 
subspace of C^**" where it is possible to correct all errors 
affecting at most t — [{d— 1)/2J qubits, i.e., correspond- 
ing to those Kraus operators having at most t operators 
■^ij U = 1, . • . ,?T-) different from identity. Furthermore, 
one denotes by Pc the projector onto the quantum code 
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C, and by £ a subset of the error operators {Ai}. The 
quantum code C is able to correct all errors in £ if and 
only if there exists a hermitian matrix 7 such that 



PcAIA^Pc 



(43) 



for any pair of error operators A\,Am € <£ (Knill 2002 1. 



The pair (C,€), consisting of a quantum code C and a 
vector space of error operators is called degenerate if 
the matrix 7 in (43) is singular; otherwise, (C, 2:) is said 
to be nondegenerate. By exploiting the symmetry of the 
map under isometry transformation of the Kraus oper- 
ators (see Sec. |ll|, it follows that in the nondegenerate 
case different error operators map the code in orthogo- 
nal subspaces, while for degenerate codes it may happen 
that distinct error operators transform the code into non- 
orthogonal subspaces. 

The action of the correctable error operators (£ on vec- 
tors spanning C (codewords) allows to define the CPTP 
recovery map o The composition of the re- 

(n) 

covery map <i>]j with the quantum channel gives the re- 
covered channel, 



1 k 



(44) 

The action of the recovery map may be interpreted 
as a unitary transformation involving an ancillary sys- 
tem, followed by a measurement on such ancilla aimed 
to extract the error syndrome and by a subsequent error 
correction action on the system. 

An upper bound on the rates achievable by non- 
degenerate quantum codes is known to be the quantum 



version of the Hamming bound ( Ekert and Macchiavello 
19961) 



2'=^3* 

■i=0 



< 2" 



(45) 



In the framework of independent and identical dis- 
tributed (i.i.d.) errors, no codes are known to violate 
such a bound. There are also upper bounds that ap- 
ply to all quantum codes, not just non-degen erate ones, 
like quantum Singleton bound n > At + k (Knill and 
Laflammel 119971). 



Unfortunately the explicit construction of practical 
codes is not an easy task. Historically the first attempts 
were made following classical linear codes (Peterson and 



Weldon 1972). A great advantage of linear codes over 



general error correcting codes is their compact specifi- 
cation. Then, along this line, quasi classical [or CSS 
(Calderbank, Shor, Steane)] codes have been proposed 
(Calderbank and Shorl 19961 Steane 1996). Another way 



to construct quantum codes is to exploit the group struc- 
ture o f the set of errors £, as it is done for stabilizer 
codes (Gottesman 1997). Also this kind of codes have 



a compact description which is mirrored into an efficient 



encoding and decoding procedure. An example is pro- 



vid ed by the [[5,1,3]] code introduced in ( Laflamme et 
al. 2009| ) saturating the quantum Hamming bound ( 45 1 . 

The above methods of constructing codes can be sum- 
marized as follow. Consider a set € of 2""^^ commut- 
ing error operators of the kind — cri"''*cri^'\ (i — 
1, . . . , 2""'^), with a,;, Pi € {0, 1}" and with a'x,y,z repre- 
senting the Pauli operators. The set of vectors stabilized 
by the error operators forms the quantum code, namely 



x) e 



ei\x) = \x),\/i = 1, 



}. (46) 



The error operators can be represented by 2n dimen- 
sional vectors v^- = {ai\/3i) in the field Fj". Let us take 
n — k such vectors named g^, . . ■g"'~^ that are linearly 
independent and write down the following (n — k) x 2n 
(parity check) matrix 
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H 



(47) 



9 



The subspace C is a [[n,fc,£i]] code, where d is the mini- 
mum Hamming distance between the vectors g^ gener- 
ating the code. The quantity Hv^, represents the 
error syndrome synd{ei). Clearly only errors such 
that synd{ei) 7^ can be identified and the set € is cor- 
rectable by C if and only if synd{ei) ^ synd(ej) for all 
i,i = 1, . . . , 2"^*^ (which is the condition corresponding 



to (|43|). 

For such codes there exists a lower bound on the rates 
(the quantum version of the Gilbert- Varshamov b ound) 
telling us that good codes exists (Calderbank et al. 1997) 



2t . . 

4=0 



(48) 



To study the performance of a given code (even in the 
presence of correlated errors), it seems natural from (41) 



to use the fidelity \A3,\ as an indicator of its effectiveness. 
However, following I Nielsen] |1996[ one is led to consider 
entanglement fidelity (see Appendix |A]) as the most ap- 
propriate quantifier of the effectiveness of quantum error 
correct ing codes (especially in the prese nce of correlated 
errors ( Cafaro et al.\ 2011)). From Eq. (A7) one obtains 
the entanglement fidelity associated to the error correc- 
tion code: 



1 

22^ 



E |Tr(i?i^k) 



(49) 



l,k 



1. Performance achieved in the case of correlated errors 

Although the idea of i.i.d. errors was underlying the 
standard theory of quantum error correcting codes, re- 
cent studies try to characterize the effect of correlated 
noise on the performance of most relevant QECCs. 
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To get some insights let us start by considering a mem- 
ory channel described by the Markovian map ( 30 1 with 
.)AI being the Kraus operators Ai^ 



unitary and 
so to write 



(1 - ^J.)p^^ +fld,, 



1, 



{A,^(»...(E)A,Jp{A, 



The correlation parameter /i describes the degree of mem- 
ory of the considered channel. For = one obtaines the 
case of i.i.d. (memoryless) noise, while the limit /i = 1 
describes completely correlated errors. 

In Ref. (Cafaro and Mancini 
Ao 



VZ, Po 



2010a|) the case with 
p and pi — p 



1 



(To = I, Ai = 
has been considered. Since this situation is analogous to 
the one with classical correlated errors, a [n, 1,3] repeti- 
tion code (with n > 3) has been analyzed. It has turned 
out that the entanglement fidelity (49 1, while obviously 
increasing with n, is quickly decreasing as fi becomes 
greater than zero. The analysis was subsequently gener- 
alized (Cafaro and Mancinil 2010b) to include depolar- 



izing quantum channels with correlated and asymmetric 
errors, i.e. Aq = aa, Ai = a^, A2 ^ ay, A3 ^ a^, 
and Po, Pi 7^ P2 Ps- In this case quantum error correc- 
tion has been performed via the 



[5, 1, 3]] and [[7, 1, 3]] 
stabilizer codes (iGottesman 120091) showing that the per- 



formance of both codes is lowered by the presence of cor- 
relations. Furthermore, it has been uncovered that the 
asymmetry in the error probabilities does not affect the 
performance of the five-qubit code while it does affect 
the performance of the seven-qubit code which results 
less effective when considering correlated and symmetric 
depolarizing errors, but more effective for correlated and 
asymmetric errors. 

Another interesting model is given by the convex com- 
bination of uncorrelated and completely correlated quan- 
tum channels (see Sec. |IV.E ) 



$ (p) = (1 - A') I'^.^o (p) + ip) > (51) 



where ^f_i=o and are given by (50 1 in the limiting 



cases of /i = and p — I (corresponding to uncorrelated 
and completely correlated errors respectively). For this 
channel, with Aq = ao, Ai = a^, Po = 1 — p and pi = p, 
the performance of a [n, 1, 3] repetition code (with n > 3 ) 
has been analyzed in Ref. ( Cafaro and Mancini[ 20111. 
Also in this case the entanglement fidelity (|49[), is quickly 
decreasing by increasing the degree of memory p. 

From the above results it seems that correlations have 
always a negative effects on the standard QECC. How- 
ever it is worth considering other situations as well. For 
instance, correlated errors can arise when qubits interact 
with a common bath. In Ref. ( [Klesse and Frank 2005) 
a quasi classical code has been analyzed, with variable 
length n (number of physical qubits) and size k (number 
of logical qubits) using a spin-boson model consisting of 



n spins describing an n-qubit register coupled to a com- 
mon bosonic bath. The possibility of exchanging bath 
bosons between qubits gives rise to spatial and temporal 
correlations in the noise. The amount of such correlations 
can be controlled by the inter-spin distance. Also in this 
case quantum error correction is substantially hampered 
by the kind of noise-correlations captured in the model. 
However, one could also consider the same code when 
qubits all interact with a common spin bath, rather than 



a boson bath. In Ref. (Shabani 2008) such a situation 



has been studied using a dephasing bath interaction and 
it has been observed that better performance can be ob- 
tained in the presence of correlated errors, depending on 
the timing of the recovery. 

As one has already remarked, generally in the stan- 
dard QECC theory single-qubit-errors where only one 
qubit has undergone interaction with the environment 
are assumed to be the most common ones. More pre- 
cisely, it is assumed that the probability of k (integer 
fc > 1) errors is of order e*^, which is much smaller than 
e, given that the probability e of a single error is small 
enough. This is the actual meaning of the error inde- 
pendence condition. However, it should be noted that 
the condition of independent errors is not equivalent to 
the condition of independent quantum channels, in which 
case each qubit interacts with its own environment that 
do not interact among themselves. Although the inde- 
pendence of the qubit-environment interactions ensures 
the independence condition, the converse is not guaran- 
teed. In Ref. ( Hwang, Ahn and Hwang |2001 1 it is shown 
that even if qubits do not interact independently with 
environments, the generated errors satisfy the indepen- 
dence condition up to second order, provided that quan- 
tum bits do not directly interact with each other. Gener- 
ally, this no-qubits-interaction condition is assumed ex- 
cept for the case where two-qubit gate operation is being 
performed. Thus, in such cases standard QECCs work 
enough well. A pictorial representation of these differ- 
ent kinds of system-environment interation is depicted in 
Fig.i 

Let us see why a QECC that corrects one qubit error 
in the standard way may fail in correcting strongly corre- 
lated errors like controlled bit flip errors. Assume to have 
a quantum code C, with orthogonal codewords {|0i)}iJi, 
that can correct any error operator belonging to € when 
applied to any vector S C. Specifically, the code 
works in such a way that for any noise operator A G € , 
a corrupted codeword = is mapped to a prod- 

uct state (S) \synd{A)) , where synd(A) stands for the 
error-syndrome of A. Since synd{A) does not depend on 
\(f>i) , one can correct errors applied to any linear combi- 

nation J2i^i\'t'i) of basis vectors {l^fiilliJi. However, 
this is no longer true for controlled bit flip errors, where 
the error may depend on the specific codeword 
In that case the corrupted codeword reads 
Then, upon decoding one gets ^iCi\(j)i) (S) \synd{Ai)) 
which, by tracing out the syndrome register, gives a 



15 




Environment 
Qubit 
Coupling 



4 4 6 6 





FIG. 4 A pictorial representation of different kinds of system- 
environment interactions. 



state different of the original one J^i'^il'Pi)- Actually, 
for controlled bit-flip errors, it is like the environment 
gets information about the codeword and thus corrupt 
it. This kind of errors is particularly relevant for mod- 
eling adversarial noise, where an adversary is allowed to 
decide which error operator to apply based on the specific 



codew ord it acts upon. In Ref . ( Ben- Aroya and Ta-Shma 
|2009) it has shown that no QECC can perfectly correct 
controlled bit-flip errors, although a QECC of arbitrarily 
high dimension can approximately correct them (i.e., it 
can give back a codeword 'close' to the original one once 
an error acted on it). 



B. Decoherence free subspaces 

A suitable strategy to deal with completely correlated 
errors is represented by Noiseless Codes, also known as 
Decoherence Free Subspaces (DFS). This is a passive 
quantum error correction method where the key idea is 
that of avoiding decoherence by encoding quantum in- 
formation into special subspaces that are protected from 
the interaction with the environment by virtue of some 
specific dynamical symmetry, see e.g. |Lidar and Whaley} 

mm 

To establish a direct link between QECCs and DFSs, 
one writes the evolution of the system Q density operator 
as 

(Pq) - Ttb [U (pq ®pb)U^^J2 ^^PQ^i ' (52) 



where U — e^''^-^'^+^^+^'"')* is the unitary evolution 
operator for the system-bath closed system and Ak the 
Kraus operators. 

It turns out that a subspace Hufs of "H is a DFS if 
and only if all Kraus operators, when restricted to "HdfSj 
are equal, up to a multiplicative constant, to a unitary 



transformation U, 



(DFS) 



the matrix representation of the Kraus operators is given 
by 



V 




Mk 



(53) 



where A/k is an arbitrary matrix that acts on H^pg (with 
H — "HdfsJB T~l-BFs) cause decoherence there. 



Equation (53) implies 



7ki-fc 







M^Mi 



(54) 



where 7ki = 5k5i- Applying condition (43) to the present 
setting, it follows that DFS can be viewed as a special 
class of QECCs, where upon restriction to the code space 
C = HdfSj all recovery operators are proportional to 
the inverse of the unitary Uq'^^'^ , 



(., 



(DFS)^ 

Q J 



(55) 



It is worth noticing that in the DFSs case the matrix 
has rank 1. Hence, a DFS is an example of maximally 
degenerated quantum error correcting code. 

To show the effectiveness of DFS for completely cor- 
related errors, let us consider a simple example ^^^^i m 
(50) with n = 2 and Aq — ctq, Ai — cr^, po — 1—p, pi = p 



so that 



^00 = Vp^o (8) ctq , Aoi = , 



Aio = 0, All = V(l -P) 



(56) 



Clearly the DFS is given by span{\01), |10)} where one 
can safely encode a qubit 



10, 



101) 



II, 



|10) . 



(57) 



By extending this argument, one can say that in the case 
of completely correlated errors it is possible to exploit the 
invariance of a subspace to safely encode information. 
In Ref. (Chiribella et al. 2010), it has been provided 



a generalized quantum Hamming bound for nondegen- 
erate codes, which depends on the rank of the Choi- 
Jamiolkowski operator representing the noise process and 
holds for any kind of noise. The original Hamming bound 
(45), which was formulated for the case of independent 



In this case, in a suitable basis. 



noise on the encoding systems is then recovered as a 
particular case. For completely correlated noise it has 
been shown how to exploit degeneracy to violate the gen- 
eralized quantum Hamming bound and achieve perfect 
quantum error correction with fewer resources than those 
needed for non-degenerate codes. As an example consider 
the following channel 

Hp) ^pp+ E (p^^'j- ^^^^^^ P ""^^^^^ 

i — l,....n,j>i 

+ VY,, of of p of of a«a^) p , 

(58) 
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where the input state is left unchanged with probability 
p = 1- J2z=i,...,n,3>iiPx,ij +PY,zj +PZAj), whilc it Un- 
dergoes Pauli operators a^, (Jy and on qubits i and j 
with probabilities px,ij, Px,ij and pz.ij respectively. By 
evaluating the rank of the map ( 58 1 one can see that the 



logical qubits) . Also in this case the entanglement fidelity 
(49 1, while obviously increasing with n, is increasing with 



generalized quantum Hamming bound of Ref. ( Chiribella 



et al. 2010 ) reads in this case 



/i, achieving the maximum allowed value of 1 when /i = 1. 

These results are opposite to those reviewed in the 
previous Section about repetition codes. Then, Refs. 



< 2'' 



(59) 



Then, by considering for instance k = 1 one gets n = 7 as 
smallest integer satisfying the bound. However, one can 
also construct codes with lower values for n. As matter 
of fact consider the following code 



(Cafaro and Mancini 2011) and (Cafaro and Mancini 
201ip pointed out that there is a certain threshold value 
JJFJj})) of fj, depending on p, below which the repetition 
code works better than the noiseless one and above which 
the noiseless code works better than the repetition one. 
These results are summarized in Fig. [5] 

Similar results were obtained in Ref. ( [D'Arrigo et al. 



\0l) = 1000), lU) = |111). 



(60) 



Now notice that these codewords are not affected by the 
action of az on any pair of qubits. Consequently the 
action of ax on a pair of qubits is identical to the ac- 
tion of CT.„ on the same pair of qubits. In other words of al'' , are transmitted without errors. Of course super 



2008), where it has been considered a Markovian cor- 
related dephasing channel. In the system-environment 
Hamiltonian, the l-th carrier (qubit) interacts with the 

environment by means of a single Pauli operator, ai^'K 
Dephasing channels are thus characterized by the prop- 
erty that, when N qubits are sent through the chan- 
nel, the states of a preferential orthonormal basis = 
\ji,---,jN),ji,---,jN = 0,1}, with eigenvectors 
(0 



the code is degenerate. Therefore, one has only to cor- 
rect errors due to a^ operators. This can be realized 
through a projective measurement onto the subspaces 
Soo = span{|000), |111)}, Sqi = span{|100), |011)}, 
5io = span{|010), |101)} and Sn = spanjjoOl), jllO)}. 
If the measurement outcome is 00, no errors have affected 
the qubits; on the contrary, if the measurement outcome 
is 01, errors have affected qubits 2 and 3 and can be cor- 
rected by applying there a^ ; if the measurement outcome 
is 10, errors have affected qubits 1 and 3 and can be cor- 
rected by applying there a^; finally if the measurement 
outcome is 11, errors have affected qubits 1 and 2 and can 
be corrected by applying there ax- It is clear that this 



positions of basis states may decohere, thus corrupting 
the transmission of quantum information. This model 
has been used to compare the performance of the stan- 
dard three-qubit repetition code with a two-qubit noise- 
less code. 

More generally, one may argue that, for the memory 
channel models (50), (51), there must be a threshold 



code violates the quantum Hamming bound ( 59 ) thanks 
to the invariance of the coding subspace under the action 
of pair of az which allows for perfect error correction. 



1. Threshold values for the degree of correlation 

Having seen that DFSs are suitable to encode infor- 
mation in the presence of completely correlated errors, it 
is natural to expect that their performance decrease by 
reducing the degree of errors' correlation (Demkowicz 



value fJ,*{p) that allows us to select the best code between 
the standard and the noiseless ones. 



C. Designing quantum codes for correlated errors 

The specific features of error models can be used to de- 
sign new quantum codes that better cope with correlated 
errors. 



For instance, in ( Novais and Baranger 2006 ) it is stud 



Dobrzanski, Kolenderski and Banaszek 


In Ref. ( 


Cafaro and Mancini 


2010a 


)t 



2007) 



le model of ( 50 ) 



with Aq = ao, Al = ax, Po = ^ — P and Pi — p has 
been considered using a noiseless code of variable length 
(number of physical qubits) and fixed size (number of 
logical qubits). It has turned out that the entanglement 
fidelity ( [49| , while obviously increasing with n, is also 
with fi and achieves the maximum allowed value of 1 for 
^ = 1. 



In Ref. (Cafaro and Mancini 2011) the model of (51) 
with Aq = ao, Ai = ax, Po = ^ — P and pi — p has 
been considered using a noiseless code of variable length 
(number of physical qubits) and fixed size (number of 



led the decoherence of a quantum computer in an envi- 
ronment which is inherently correlated in time and space. 
They first derive the nonunitary time evolution of the 
computer and environment in the presence of a stabilizer 
error correcting code, providing a general way to quantify 
decoherence. The general theory is then applied to the 
spin-boson model. The results demonstrate that effects 
of long-range correlations can be systematically reduced 
by small changes in the error correction codes. 

More generally, one can summarize the following de- 
velopments. 



1. Concatenated codes 

The results at the end of the previous Subsec- 
tion suggest that it may be convenient to concatenate 
decoherence-free subspaces with standard quantum error 
correcting codes in order to achieve higher entanglement 
fidelity values in both low and high correlations regimes. 
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This kind of concatenation was first introduced in 



Ref. (Lidar, Bacon and Whaley 1999), however in Ref. 



( Clemens, Siddiqui and Gea-Banacloche 
sequently in Ref. (ICafaro and Mancini 



j2004D and sub- 
20TTf it 



was in- 



vestigated in the context of memory channels. 

For the sake of simplicity, let us illustrate the case 
of two layers of concatenation and consider single qubit 
encoding. Assume the inner code (first layer) is a 
[[ni, /ci, di]] stabilizer code Ci with generators Gi = 
{gj : i = 1,..., ni — ki} , and the outer code (second 
layer) is a [[n2j Ij c?2]] stabilizer code C2 with generators 
G2 = {9j ■ j = ^T--, nj — 1]. The concatenated code 
C — Ci o C2 maps ki qubits into n — nin2 qubits, with 
code construction parsing the n qubits into ni blocks 
B (b) (b = 1,..., rii) each containing 712 qubits. In other 
words, given a codeword |cin) for the inner code Ci, 



|Cin) = I'^j) ' 

J=0 



(61) 



where {|0j)} are basis vectors for Ci, the concatenated 
code C is constructed as follows. For any codeword |cout) 
for the outer code C2, 

|cout) = ^ ail,..., i„J«i---»n2) , (62) 
n,..., i„„ 



with — ®.. 

ki — 1 for j = 1 
in Ci, that is 



with = 
vector \(, 



„2 ) , replace each basis vector 
n2 by a basis 



|0,,)(8)...«)|<^i„J. (63) 



Further details on the construction of the stabilizer gen- 
erators of C can be found in (Gaitan 2008). As a fi- 



nal remark, one simply points out that the above men- 
tioned construction produces a [[niri2, fci, d]] code with 
d> did2- 

As an illustrative example, consider to encode one logi- 
cal qubit with a concatenated subspace obtained by com- 
bining the decoherence free subspace (for a completely 
correlated bit flips) 



|0l) = 



II/,) = 1- 



(64) 



(inner code, Cdfs — Cinncr ) with the three-qubit bit 
repetition code in ([6O]) (outer code, Cbit — Coutor)- One 
obtains that the codewords of the concatenated code C = 
Cdfs ° Cbit are given by, 



|0l) 



(1000000) - 1000111) + 1111000) - 1111111)) : 
(1000000) + 1000111) - 1111000) - liiiiii)) . 



The entanglement fidelity for the concatenation of a 
repetition code and a noiseless code for the model of ( 50 ) 
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FIG. 5 Entanglement fidelity for the model (50 1 top figure 
and (51 1 bottom figure witli p = 10~^. Blue solid lines corre- 



spond to standard code. Red solid lines correspond to noise- 
less code. Dashed lines correspond to their concatenation. 



with ^0 = Co, A\ = ax, Po = ^ — P and pi = p are 
reported in Fig. [5j It turns out that the concatenated 
code does not work well for partially correlated errors. It 
is always better to use either the outer or the inner code 
alone depending on whether one is below or above the 
threshold value /i*(p). 

The entanglement fidelity for the concatenation of rep- 
etition code and noiseless code for the model of ( 51 ) with 



CTqi Ai = (Tx, Po — 1 ~ P and pi — p are shown in 
Fig. [5] too. Here the concatenated code works perfectly 
almost everywhere. Hence, one may argue that for the 



model of (51 1 the concatenation trick results particularly 
advantageous in the presence of partially correlated er- 
rors. 



2. Burst errors codes 

In communication channels noise disturbances over 
long time periods lead to temporally continuous errors, 
as opposed to noise disturbances over short time periods 
leading to random single errors. Analogously, stored data 
can be defective over spatially length scales greater than 
a single bit. Such errors are called burst errors. They 
are well studied in the classical framework where cor- 
responding error correc t ing c odes h ave been developed 
(Peterson and Weldon 1972). In (Vatan, Roychowd- 
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hury and Anantram 1997), quantum analog of burst- 



error correcting codes have been considered. Hamming 
and Gilbert- Varshamov type bounds have been derived 
showing that these codes are more efficient than codes 
protecting against random errors. In fact, to protect 
against burst errors of width b (that is, errors occurring 
on a number b of consecutive qubits with b a fixed con- 
stant), it is enough to map n — logn — 0{b) qubits to 
n qubits, while in the case of t random errors at least 
n— t log n qubits should be mapped to n qubits. Based on 
binary cyclic codes explicit constructions have been pre- 
sented of these almost optimal quantum code for bursts 
of width not bigger than 4 and for specific values of n. 
It remains open the problem of how to generalize these 
constructions. 

Quite generally, one must increase the redundancy of 
the codeword with increasing length of the burst. On 
the other hand, if details of the decoherence mechanism 
of qubits arc known, then it might be possible to design 
a more efficient error-correction method than the con- 
ventional scheme. In (iKawataba 2000), it is proposed a 



method for constructing quantum burst-error-correcting 
codes from known quantum error-correcting codes. This 
method is based on the interleaving technique. By using 
this method, the quantum codewords can be distributed 
amongst the qubit stream so that consecutive words are 
never next to each other. On de-interleaving they are re- 
turned to their original positions so that any errors that 
have occurred become widespread. This ensures that any 
burst (long) errors now appear as random (short) errors. 



3. Convolutional codes 



It is possible ( Chau 1998 1999 OUivier and Tillich 



2003 2004) to extend the notion of stabilizer codes to 
codes which allow for an overlap between the individual 
steps of the encoding operation. It can be noticed that 
their encoding op erations can be d escribed as quantum 
memory channels (Giitschow 2010), due to the fact that 



some of the output qubits of the nth encoding step will 
be used as inputs in the (n -I- l)th step of the encod- 
ing, thus the blocks overlap. These qubits correspond 
to the memory system of the memory channel describing 
the encoding. The encoding operation is the same in ev- 
ery step (neglecting initialization and finalization), thus 
every step is described by the same channel. From the 
classical coding theory, codes with these properties are 
called convolutional codes. 

In particular, convolutional stabilizer codes correspond 



to Clifford memory channels (Giitschow et al. 2010). 



Looking at the encoding operation as a quantum chan- 
nel, it is convenient to take the encoding of one block 
as one use of the channel. As one passes on qubits from 
one step of the encoding to the next, one deals with a 
memory channel. All encoding steps are the same, so 
our channel can be concatenated and used for arbitrary 
many steps of the encoding procedure. Stated more pre- 



cisely, for every [[n, fc,m]] convolutional stabilizer code 
one can find an encoding operation which is described 
by a the concatenation of a reversible Clifford memory 
channel. The channel has n input and output qubits and 
uses m qubits of memory. In this framework the memory 
system appears as a resource for the implementation of 
the correcting code. 

Unfortunately convolutional codes also carry disadvan- 
tages. Because information is transmitted from one block 
to the next, errors can spread as well. Depending on the 
encoding algorithm catastrophic errors can occur. Catas- 
trophic errors happen during the process of encoding (or 
transmission or decoding) and affect only a few qubits 
at first. During the encoding (or decoding) they spread 
without bound. These errors can only be corrected with 
operations of unbounded support (i.e., they act on an 
unbounded number of qubits). If one considers trans- 
missions of finite length, the support of the operation is 
only bounded by the length of the transmission. There is 
one exception of this behavior. The support of a catas- 
trophic errors can also stay finite but the localization 
area moves in an unbounded manner, i.e., the error will 
always affect the qubits at the end of the transmission. 
In both cases the error can only be corrected at the end of 
the transmission. This implies that the decoding process 
can start only after the end of the transmission. There- 
fore the delay between receiving the first qubit and de- 
coding it is as long as the transmission itself. Further- 
more the memory requirement on the decoder side is of 
the order of the length of the encoded message. Con- 
sequently these errors have to be avoi ded by employing 
non-catastrophic convolutional codes (Grassl and Roet 



teler 2006). Non-catastrophic codes with minimal di- 



mension of the memory system can be devised as well 
(Houshmand, Hosseini-Khayat and Wilde 2011 Wilde 



et al. 2011). Moreover, the dimension of the memory 



system of non-catastrophic codes can be further opti- 
mized by allowing the assistance of pre-shared entangle- 
ment (|Wilde and Brunl |2009l |2010a|b[ ). The link be- 



tween Clifford memory channels and convolutional stabi- 
lizer codes allows to prove that a convolutional encoder 
is non-catastrophic if and only if the corresponding mem- 



ory channel representing is strictly forgetful (Giitschow 
20101). 



Not only convolutional codes can be described by mem- 
ory channels, but also memory channels could benefit 
from convolutional codes. In fact, when sending informa- 
tion across correlated channels, data put into one chan- 
nel may influence the output of all further invocations. 
Therefore, a decoding circuit needs to gather information 
from the output of several channel uses in order to de- 
code any given bit. Block codes are ill- adapted to this 
task, as one expects a considerable loss of information 
to occur at every block boundary. Convolutional codes 
could be much better suited for this particular situation. 
Lacking the boundaries between code blocks, convolu- 
tional codes exhibit the same "continuous structure" as 
channels with memories. Convolutional decoding circuits 
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continuously gather output from the channel, outputting 
decoded qubits as soon as sufhcient data has become 
available. 



VI. CODING THEOREMS AND CAPACITIES 

Here, one discusses more in detail and quantify the ca- 
pability of a quantum channel $ to transfer classical and 
quantum information encoded on quantum systems. This 
will lead to the notions of capacities of a quantum chan- 
nel, which are the quantum counterparts of the Shannon 



capacity (Shannon 1948) 



In the context of classical communication theory, the 
latter was defined as the optimal rate of reliable infor- 
mation transmission, that is the maximum rate at which 
information can be sent through a channel with asymp- 
totically vanishing error probability, in the limit of long 
messages, and by using proper encoding and decoding 
procedures (asymptotically eliminating all errors) . Shan- 
non proved that this quantity is expressed in terms of a 
maximization of an entropic quantity (mutual informa- 
tion) between the two parties. 

Being non constructive, this proof does not provide an 
explicit recipe for the optimal code. Roughly speaking, 
the codewords of the optimal code correspond to distin- 
guishable output states, in such a way one can efficiently 
recover the input with an arbitrarily small error. 

The generalization of this concept to the quantum set- 
ting leads to the definition a family of channel capacities, 
depending on whether classical or quantum information 
has to be transmitted, and whether additional resources, 
as pre-shared entanglement, are exploited. Indeed, the 
capacities of a quantum channel are operationally de- 
fined as the optimal communication rates through the 
given quantum channel. Then, coding theorems prove 
their expression in terms of entropic quantities (see Ap- 
pendix Ib|) . 



A. Operational definitions and memoryless setting 

The classical (quantum) capacity of a quantum chan- 
nel is then the maximum rate at which classical (quan- 
tum) information, encoded on a quantum system, can 
be reliably sent from the sender to the receiver through 
the given quantum channel. This rate basically is the 
maximum ratio among the number of bits (qubits) trans- 
mitted and the "redundancy" employed into the code. 
In mathematical terms, the classical capacity C($) of a 



noisy quantum channel can be defined as follows ( Ben- 
nett and Shorl 119981) 



C($) := linilimsup |- : \$(^^-..)^^(^^^.) V|^)g{|o),|i)}«. 
such that 

F(|V')(V'|,<I'^'"'^°<i>®"o4^-")(|^)(V.|)) >l-e} , 



with F{pi,p2) being the quantum fidelity between the 
states pi and p2 as defined in Eq. (A3) in Appendix |A 
$ E being an encoding map from n qubits to m inputs o ' 
the channel $, and $_d a decoding map from m channel 
outputs and n qubits. Let us point out that the input 
state IV') is simply a tensor product of states |0) and 
|1), hence encoding classical information. Indeed, the 
definition of C($) is related to the transmission of clas- 
sical information encoded in quantum states, hence no 
superposition states are considered here. In other words, 
C{^) is defined as the optimal rate (i.e., transmitted bits 
per channel use) at which the sender can send a tensor 
product state, i.e. \tp) G {|0), ll)}*^", of n qubits, for ar- 
bitrarily large n, to the receiver, who is able to recover it 
with fidelity greater than 1 — e, with arbitrarily small e, 
after encoding, transmission, and decoding procedures. 

Another quantity of interest is the one-shot capac- 
ity, or product- state capacity, denoted Ci($), which 
is defined as in ( 65 ) with additional assumption that 



the employed codes are made of separable states, that 
is, <i>^~*'"^(|V')(V'l) is a separable state for all £ 

{|o),|i)n. 

Similarly, the quantum capacity (5($) is defined as 
(IBennett and Shorl 119981) 



g($) linilimsup - : 3^ <j(fc^") ^(-^fc) V|^)g^»fc 

<^^^ fe— >oo ' ^ ' ^ 

such that 

F (|V')(^|,<f^"^'' o o $(,^^")(|^)(^|)) > 1 - e} , 

(66) 

where Ji is the two-dimensional Hilbert space of a two- 
level quantum system. Notice that the only difference 
between the definitions of C and Q is that now the in- 
formation is quantum and can be encoded also on su- 
perpositions of k qubits. Since a tensor product state is 
always a specific case of a state living in Ji, one obtains 
the first of a series of inequalities between channel capac- 
ities, i.e. Q{^) < C(^). This crucial difference between 
C and Q allows one to exploit some peculiar properties 
of quantum states, as entanglement, in order to build 
up optimal encoding inputs. Analogously, joint quantum 
measurements of the output states have to be taken into 
account in the choice of the optimal decoding strategies. 

The entanglement-assisted classical capacity, denoted 
Cea($), is defined as the maximum rate of reliable trans- 
mission of classical information when the sender and the 
receiver have a their disposal an unbounded number of 
pre-shared maximally entangled states. Coding theorems 
prove that the entanglement-assisted classical capacity of 
a memoryless channel is given in terms of the quantum 
mutual informa tion [d efined in Eq. ( B4 1 in Appendix [b] 
(IBennett et aLl 119991 120021): 



Cea{'i>) = maxl(p, $) 
p 



(67) 



(65) where the maximization is over any input p. 
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Concerning the (unassisted) classical capacity, let us 
first focus on the one-shot classical capacity, whose ex- 
pression is given by ( [Holevo , 1998 Schumacher and West 



moreland 1997) 



Ci{^) = max 
{pi,Pi} 

i t _ 

(68) 

where the maximum is taken over all probability distribu- 
tions {pi} and collections of density operators {pi} (pos- 
sibly satisfying some additional input constraints). The 
quantity in bracket in Eq. ( 68 ) is known as Holevo in- 



formation and denoted as x(<I>, {pi, Pi}). The expression 
for the classical capacity of the quantum channel is then 
ob tained by allowing encoding blocks of increasing length 
n (IHolevol 119981 ISchumacher and Westmorelandl 119971), 



I.e., 



C($) = lim -Ci(4>®"). 



(69) 



Since an additional resource (shared entanglement) is 
required for the definition of Cea{^), the following in- 
equality holds C($) < Cea(<I>). Finally, let us point out 
that the additivity conjecture, C($) — Ci{^), i.e. en- 
tangled input states do not increase the classical capac- 
ity, was believed to be true, allowing a practical compu- 
tation of this quantity. However, recently, it has been 
shown that this conjecture is not always true, by con- 
structing a counter-example where the one-shot capacity 
Ci is superadditive, i.e. Ci(<I> (g) $) > 2Ci(<I>) with 



$ being a random- unitary quantum channel (Hastings 



2009 ). Notwithstanding, one may safely say that additiv- 



ity holds true for a large variety of channels ( Amosov and 
Mancini[ [2009| [Hiroshima! |2006[ ). Note also that the ad- 
ditivity counter-example does not necessarily imply that 
C($i (g) <I>2) > C($i) + C($2) since an additional regu- 
larization is needed in computing C as in Eq. (69). These 



results leave still open the problem of channel capacity, 
i.e. the full optimization over all entangled input states 
is generally required to calculate the maximum capacity 
C. 

The expression for the quantum ca pacity Q(^) of a 
memoryless channel $ is provided by (Barnum, Nielsen 



and Schumacherl [19981 fPevetakt [20051 |Lloydl |1997[ |Shor 
200l 



g($) = lim -Qi($) 

n— >-oo 72 



(») 



where 



Qi($)(") =maxJ(p,$^") 



(70) 



(71) 



where the maximization is performed over all input states 
for n uses of the quantum cannel and J{p, $®") is the 
coherent information of the map <I)®", as defined in 
Appendix [Bj The maximum of J may be strictly su- 
peradditive for parallel channels, hence the quantum 



capacity may be greater than its one-shot expression, 
i.e. Q($) > Qi($), in contrast to the case of the 
entanglement-assisted classical capacity Cea which is al- 
ways additive. For this reason, the regularization for n 
which tends to infinite (i.e., lim„_>.oo) is necessary. On 
top of that, the coherent information is not, in general, 
a concave function, therefore there might be local max- 
ima which are not global ones and this makes the cal- 
culation of Q one of the most difficult task in quantum 
communication theory. Notice that for degradable chan- 
nels, defined in Sec. [Hi the coherent information can be 
proved to be additive and then the regularization for n 
is not necessary, i.e. Q{^) = Qi{^) — maxp J{p, $). On 
the other side, for anti-degradable channels $ one always 
has (5($) = 0. The latter holds also for the class of PPT 



maps (see Section II. A). The three channel capacities. 



defined so far, satisfy the following inequality 

Q(*) < C($) < Ceai^) . 



(72) 



Finally, recently an example of a new interesting phe- 
nomenon, called superactivation, has been provided 
( Smith and Yard[ 2010[ ). In particular, it has been shown 
that it is possible to find channels $1 and $2 with zero 
capacity, i.e. C($i) = Q($i) = C($i) = Q($i) 0, for 
which, by joining the two communication lines $1 and 
$2 in $1 (81 $2, it becomes possible to transmit quantum 
information, i.e. Q($i $2) > 0. Specifically, $1 and 
$1 are given by an anti-degradable channel and a PPT 
one, - see ([Brandao, Oppenheim and Strelchuk[[2012[) for 



other examples in terms of depolarizing maps and for a 
more general construction. 



B. Entropic upper bounds: memory setting 

The operational definitions of channel capacities, in- 
troduced above for the memoryless setting, apply also 
straightforwardly to memory quantum channels. Indeed, 
they are quite generally related to the optimal classical 
and quantum information transmission rates between two 
parties, no matter how complex is the internal structure 
of the channel. However, the memory setting is often 
more complicate than the memoryless one since addi- 
tional channel features have to be specified in order to 
discuss the concept of communication capacity. For in- 
stance, in the case of non-anticipatory channels in Sec. 
III.B.1[ the map can be described in terms of an addi- 
tional memory system Af , and then the initial and final 
memory states have to be characterized before defining 
the channel capacities. In particular, one has to distin- 
guish different setups where Alice, Bob and/or Eve (third 
party) may control (or use for the encoding/decoding 
procedures) the initial/final states of the memory system 
( [Kretschmann and Werner[ 2004). Hence, this leads to 
different definitions and bounds for the memory channel 
capacities. In fact, for a channel $ one can at least de- 
fine four classical capacities C^b($), Cae{^), CEB.fii^): 
CEE,fj,i^), with the first (second) index representing the 
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party controlling the initial (final) memory state and fi 
being the Eve's choice for the initial state of M (when 
considered). The same classification holds for the quan- 
tum capacities Qab{^), Qae{^), QEB,ti.{^), Qee,,i{^)- 
Moreover, it may be generalized to setups where addi- 
tion resources, as shared entanglement and classical side 
(forward and/or backward) communication, do assist the 
memory noisy transmission. 



1. Finite-memory channels 

The first attempts to extend the notion of classical and 
quantum capacity to memory channels were shown in 
Ref. (Bowen and Mancini 2004), where some bounds in 



terms of mutual and coherent information were provided, 
independently from the internal evolution of the mem- 
ory state. In particular, they considered finite-memory 
channels, i.e. maps with memory of finite dimension (de- 
scribed in Sec. III.B.2), exploiting the Holevo quantity 



as an upper bound for the classical capacity C(<i>), i.e. 
for a given initial memory state 

C{^) < lim - max x{^^"\ {p^ , P^}) (73) 

n— >-oo n {pi,Pi} 



where x(. . . ) is the Holevo quantity defined in Eq. (68 1, 
as in Eq. (191, and pi denotes a set of Alice states 



(shared entangled states between Alice and Bob) for 
unassisted (entanglement-assisted) communication. Un- 
likely in the memoryless case, here this upper bound is 
not always achievable. Indeed, the coding for each chan- 
nel cannot be divided into smaller blocks because the 
memory state can be entangled over several blocks. How- 
ever, this bound turns out to be achievable for a class of 
Markovian channels, i.e. maps represented by unitary 
Kraus operators and steady-state probabilities as initial 
error distributions, as defined in Sec. IV. C (Bowen and 



Mancini 2004). The intuitive explanation for this bound 



achievement is that the memory plays a small role on 
the asymptotic behaviour of the channel. In other words, 
for Markovian memory channels the capacity is only af- 
fected by the loss information to the external environ- 
ment, asymptotically neglecting its loss into the memory 
system. Analogously, one can derive a bound for the 
quantum capacity in terms of the coherent information, 
i.e. 



g($) < lim -max J(p,$(")) 



(74) 



2. Perfect memory channels 

These maps can be described in terms of Kraus op- 
erators whose number of elements is sub-exponentially 
growing in the number of channel uses n - see more de- 
tail in Sec. |IV.B[ In other terms, the size of the envi- 
ronment is not large enough to include the information 



sent from Alice to Bob, which is exponentially increas- 
ing in n, and so asymptotically the lost of information 
into the environment is negligible. This intuitively ex- 
plains that perfect memory channels are asymptotically 
noiseless and have maximal capacities. Specifically, one 
can verify that the following theorem holds (Giovannetti, 
IBurg arth and Mancini| |2009[ [Kretschmann and Werner[ 
|2005] : 



Theorem 1 L et C be a multi-use perfect memory chan- 

n g N} be the sequence 



net as m 



77|j and let {4"^ 



satisfying Eq. {29). Then for sufficiently 



exists a zero-error classical code C of size 



\C\> 



the 



(75) 



with a rate Rc > log c?q — ^ log c?^'' that for n 
converges to the optimal value log d, 



Here, dq''"'^ is the 



size of the system carriers, i. e. d, 



(«) 



dimH 



(") 
Q ■ 



Analogously, for sufficiently large n there exists a zero- 
error quantum error correcting code Q of size 



|Q|> 



(4 



(76) 



with rate Rq > logdg - i log[(4"0^ + (4"0'] which, 
again, for n ^ oo converges to the optimal value log dq . 



3. Non-anticipatory memory channels 

For generic non-anticipatory (n-fold conc atenated) 
memory quantum channels <I'^"'-' as in Eq. (17), the clas- 



sical capacities are upper bounded as ( Kretschmann and 



Werner 2005): 



C_E_E,/x($) < limsup- max x(trM o Pj) • 

n-i-oo n {pi.Pi} 

(77) 

Analogously, one finds the following upper bounds for 
the corresponding quantum capacities, i.e. 

Qee,ij.{^) < limsup — max J(trMO$("),M(»p) . (78) 



Similar expressions are obtained in the other three con- 



figurations discussed in Sec. VLB The proof of this the- 
orem comes from similar bounds for the memoryless set- 
tings. However, in this context coding theorems seem to 
be much harder to be proved, and these entropic bounds 
for classical and quantum information transmission rates 
are demonstrated to be achievable for a specific class of 
memory channels, called forgetful ones, as further dis- 
cussed below. 
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4. Forgetful channels 

These communication channels are characterized by 
the property that the effect of the initial memory state 
becomes negligible with time, i.e. memory effects die 
away exponentially fast, as discussed in Sec. |IV.D| This 
feature allows one to prove that the upper bounds in 
Sec. VLB. 3 can be actually asymptotically achieved 



(Kretschmann and Werner, 2005). This important result 



can be demonstrated reducing the problem into the case 
of memoryless channel by means of a double-blocking 
procedure. In other words, one can define the memory 
depth of a forgetful channel $ as the smallest integer 
m £ N such that 



^ r, 



(7'(g)idg™)o$ 



with V being the completely depolarizing channel, and 
idQ being the identity map acting on @{Hq). It a fi- 
nite m does exist, these maps are called strictly forget- 
ful. Hence, it is possible to group the channels into 
blocks of length m + I, encoding the input in the I chan- 
nels and ignoring the intermediate m ones. The asymp- 
totic behaviour of the channels is obtained when / — >■ cxd. 
Furthermore, restricting the inputs to product states of 
block length m -I- / and exploiting the strict forgetfulness 
property, the output state factorizes, and the whole map 
is reduced to a memoryless channel on the larger input 
space 'Hq'^'"'"™. It allows to extend Holevo and Deve- 
tak coding theorems for memoryless channels to forgetful 
memory channels. Note that the double-block strategy 
can be applied even if the channel is not strictly forgetful. 
Therefore, the memoryless expressions for classical and 
quantum channel capacities in Eqs. (69 1 and (701 can be 



applied also in the memory setting for forgetful maps, 
and the entropic upper bounds in Eqs. ( 77 1 and ( 78 1 are 
exactly achieved. 

Forgetful channels have been proven to constitute a 
dense set with th e topology induced by the norm o f com- 
plete boundness (Kretschmann and Werner 2005). That 



implies that any non-forgetful channel can be approx- 
imated by a forgetful one. Notwithstanding, their ca- 
pacities may be different. An example is given by in 
the context of Markovian-correlated channels. A long 
term-memory channel, Sec. |IV.E| , can be approximated 
by a forgetful Markovian-correlated channel. Sec. |IV.C| 
However, according to the coding theorem for long-term 
memory channel reviewed in the following section, the 
capacity of the latter does not approximate the capacity 
of the former. 



5. Long-term memory channels 

To provide coding theorems for channels with long- 
term memory (i.e., channels which are "not forgetful") is 
more complicate. A working class of this kind of quantum 
channels has been identified in ( jPatta and Dorlas[ |2007| , 
where convex combinations of memoryless channels were 



considered (see Sec. IV. E). For a channel $ in this class, 
: ©f-H®") ^ eijC^) and the action of on any 
state p*^"' e ©('H®") is given as follows 



M 

$(»)(p(")) = ^p,$f"(/9(") 

1=1 



(80) 



where <^>^ : &{H) &{K.) {i = 1,...,M) are com- 
pletely positive, trace preserving (CPT) maps andpi > 0, 

J2i=iPi — 1- Here, V. and K. denote Hilbert spaces. A 
classical version of this channel was introduced by Ja- 



cobs ( 


Jacobs 


1962 


( Ahlswede 




19681. 



The following coding theorem for quantum channels 



(79) ([80]) has been derived in patta and Dorlas^ 2007) ) 



Theorem 2 The one-shot capacity of a channel $, with 
long-term memory, defined through iSO^, is given by 



Ci($)= sup min xi<^>^AP J, Pj}) , (81) 
{pj-.Pj} ' 

The supremum is taken over all finite ensembles of states 
Pj G B{'H) with probabilities pj . 

This theorem has been proved by employing a quan- 
tum version of Feinstein's Fundamental Lemma (Fein- 



stein 


1954 


Jacobs 


19571 and a 


Strom 


s theorem ( He 


Strom 




19761 



For a quantum mem- 
oryless channel, the method yields an alternative proof of 
the Holevo-Schumacher- Westmoreland (HSW) theorem 



(Holevo 1998 Schumacher and Westmorelandl 1997) 



similar to the proof in ( [Winter 1999 ') 



The one-shot capacity can be generalized to give the 
classical capacity of the channel $ in the usual manner, 
that is, by considering inputs which are product states 
over uses of blocks of n channels, but which may be en- 
tangled across different uses within the same block. The 
classical capacity C(<I>) is obtained in the limit n — oo 
and is given by 



C($) = lim -Ci($(")). 



(82) 



These results have been generalized to the case of a con- 
tinuous set of channels by [Bjelakovic and Boche 2009[ 

Then, a coding theorem for communication assisted by 
entanglement in the c lass of channels specified by (|80t h as 
been derived as well ( Datta, Suhov and Dorlas 2008 1 . 



Theorem 3 The entanglement-assisted classical capac- 
ity of a channel <I>, with long-term memory, defined 
through (Sdj^ , is given by 



Ceai^) = max minT(p, $i) 

peeiUA) L i 



(83) 



where T{p,^i) is the quantum mutual information {B4-) 
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The proof of this theorem, makes use of the expression 



for the product state capacity (81). 

The above theorems can be extended to quantum chan- 
nels with arbitrary Markovian correlated noise as shown 



in (|Datta and Dorlas 


20091) 


and Notzel 


2009 


I provided 



sion for the the quantum capacity of long-term memory 
channels. 



6. Ergodic cq-channels with decaying input memory 

For ergodic cq-channel with decaying input memory 
a coding theorem has been derived in (Bjelakovic and 
Bochel[2008|. 



Theorem 4 Let W : x B"^ —?' C be a stationary er- 
godic decaying input memory cq-channel, then the classi- 
cal capacity is given by 



C{W)= sup iip,W), 

p stationary ergodic 

where 

W) hm - {S{p;) + S{p'^) - 

with 



(84) 

(85) 

(86) 
(87) 



Here p^^ denotes the density operator of the output state 
VF"(a;", •), x" € ^" and |a;") = e^^^ . . . (8i e^^ for some 
orthonormal basis {e^}!^']^ o/C''^'. 



The sup in Eq.(84) is calculated over all stationary 
ergodic probability measures p on A^. That is, consider a 
shift T : A^ — ^ of double infinite sequences of A, then 
p is stationary if p{Ta) = p(a) for all aA^. Moreover, it 
is ergodic if for all aA^ such that Ta = a it is p{a) = 
or 1. 

The above theorem results as an extension of cod- 
ing theorem for input memoryless cq-channel whose 
proof combines Wolfowitz's code construction ( Wol- 
fowitz[ 1957 ) and a version of the Feins tein's lemma 
( [Blackwell, Breiman and Thomasi"aii |1958 ) based on the 
notion of the joint input output probability distribution. 



VII. SOLVABLE MODELS OF MEMORY CHANNELS 

Up to date, only few models of memory quantum chan- 
nels have been fully solved in terms of capacities. One is 
the dephasing channel (in the discrete variable setting), 
and the other is the lossy bosonic channel (in the contin- 
uous variable setting) , with different types of correlations 



Oh 

•z 



LE- 
C, 

Ci 









< 


» 

( 


» 



<£r, 



H 

H 
O 



FIG. 6 Circuit representation of Eq. (91 1. The local environ- 
ment LE interacts through the unitary couplings Uj with one 
carriers at a time. Between two consecutive interactions with 
the carriers instead LE undergoes the dissipative evolution 
described by the transformations ■ 



A. Discrete memory channels 

The correlated dephasing channel was first discussed in 
( Giovannetti 2005 1 as an example of dynamical memory 
model. Such a model consists in a communication line 
where messages are encoded into identical carriers qi, q2, 
■ ■ ■ propagating through a medium E that separates the 
sender (Alice) from the receiver (Bob). 

The carriers are organized in a time-ordered sequence 
s — {ti, T2, - • •} with Tj > being the time interval be- 



tween the instants ij+i and tj at which qj+i and qj enter 
E respectively. The effective transit time is assumed to 
be constant and much shorter than the intervals tj. 

Two distinct components of the medium E are iden- 
tified: a finite dimensional Local Environment (LE) 
component which is directly coupled with the carriers 
through unitaries Uj, and a huge Reservoir (R) compo- 
nent which is coupled with LE but not with the carriers. 
The action of R on LE is supposed to induce a dissipative 
dynamics which transforms any initial states of LE into a 
stationary configuration do, with te being the character- 
istic time of the process. This is equivalent to introducing 
a one-parameter family {£t}t>o of CFTP maps which, 
given (T the initial state of LE at some time tg, repre- 
sents its evolution at time to -\-t with the density matrix 
£r(o'). In this formalism Sq coincides with identity map 
on T-Lle- On the other hand the stationary state ctq of 
LE is defined by the property 

£r{<yo) = Vr>0, (89) 

while the characteristic time te by the property 

^r>r^(0) =^oTr{e} , (90) 

for all bounded operator O of T-Lle- 

The resulting transformation in "He ® T~Lle is a time 
ordered product of interactions Uj and relaxation pro- 
cesses (see Fig. l6| . Assuming LE to be initially in the 
stationary state ctq this gives a discrete family of CFTP 
maps {"I'i"''}^ for an n carriers state p^^\ 

pin) ^ $(")(p(")) = TrLE {Un O 5.„_, oUn-1 O 

■■■oS^^oUi (p(")®fTo)},(91) 
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where I4j{- ■ • ) stands for the unitary mapping Uj{- ■ ■ )Uj 
on Ticj <8) Hle, and "o" indicates the composition of 
super-operators. 

Because of the time ordering of Eq. (91) the output 



state of a carrier might depend on the input state of the 
carriers which precedes it in s but it is always indepen- 
dent from the input state of the carriers which follows it 
in the sequence. As a matter of fact Eq. (91 1 exhibits the 



general structure of causal memory channel with time- 
translation invariance see Sec. III.B.l and (Kretschmann 



and Werner 2005) 



Now, assume Alice is producing a sequence s with 
intervals Tj greater than or equal to the characteristic 
relaxation time te of the dissipation process. In this 
case, after each interaction, the local environment LE 
has enough time to relax into the stationary configura- 
tion (To before a new carrier begins interacting with it. 
Under this hypothesis Eqs. (90) and (91) yield 



(92) 



where Af is the CPTP map which transforms the density 
matrices p of a single carrier into 



(93) 



Equation (|92| describes a memoryless configuration 
where the noise acts on the qj independently. 

Alternatively, suppose that only a subset Qm of the 
transmitted carriers is used to encode messages to Bob. 
The remaining carriers (subset Qp) are instead employed 
for perturbing LE in such a way that the Qj on which the 
messages are encoded have a better chance to reach Bob 
without being corrupted. A simple implementation of a 
noise attenuation scheme is such that the Q p carriers are 
composed by uniform strings of n states po separated by 
equal time intervals r. The information is instead en- 
coded a single carrier and the whole structure repeats af- 
ter a relaxation time te- In this configuration the trans- 
formation of the Qa carriers which comes from solving 
Eq. ( [91] ) can be computed as follows. First, one deter- 
mines the modified state (t„ of LE which arises from the 
interactions with the B carriers. This is accomplished by 
solving the set of coupled equations 



(t;. =Trc {U{po®cjj)} 

CTj + l = fr(CT') , 



J =0,1, 



1. (94) 



The density matrix cr„ which results from (94) is then 



used to determine the output state of the A carriers ac- 
cording to the equation 



Af{p)^i:TLE{U{p®an)} 



(95) 



The transformation ( 95 ) is in general different from 
Eq. (93 1 and depends explicitly on the parameters n, r 



and poj so it can be regarded as a memory channel. 

The correlated dephasing quantum channel can be con- 
sidered as an example of the above described dynamical 
model. We assume the carrier-LE interaction C/o to be 



to a control-unitary such that when the carrier is in |0)g^. 
nothing happens to LE, while when qj is in |l)q^. the en- 
vironment undergoes to the transformation 



-9 



(96) 



with g e [0, 1]. Moreover, one will assume the relaxation 
process {£t}t acting on LE to be described by amplitude 
damping maps which takes the state |1)l_e to |0)l£; with 
probability 1 — ry(r) where ?7(t) G [0, 1] is a non increasing 
function of r with characteristic time te, i.e. 

fr(|0)L£(0|) = |0)l£(0| 

£r{\l)LE{l\) = ll{T)\l)LE{l\ + {l-v{r))\Q)LEm 
£rmLE{l\) = y^|0)L£;(l| . (97) 

In this example the stationary state ctq of LE is hence 
\^)le- The parametrization of the memory effect is given 
by 77 (r), with 77 = corresponding to the memoryless case 
(fast environment relaxation, say r 3> t^) and 77 = 1 
corresponding to perfect memory case (no environment 
relaxation, say r ^ te)- Under the above conditions, it is 
possible to show that both the map M of the memoryless 
case and the map M of the noise attenuation protocol 
correspond to a phase damping channel Vg where the 
coherence terms of the input qubit p are degraded by the 
factor g — \ — 2pz with the probability of z error, i.e. 



^.(|0)c(0|) = |0)c(0|, 
^.(|0)c(l|) = 5|0)c(l| 



Va{\l)c{l\ 



\l)n{l\ 



(98) 



In particular Eq. (93 1 gives N — Vg. On the other hand, 
Eq. ( 95 ) gives M — V-g where ^ is a complicated expres- 
sion of the parameters po, and r. By appropriately 
selecting the values of the above quantities one can make 
A/" less noisy than M by having g > g. 

In the case of the phase damping channels the capac- 
ities can be explicitly computed. For instance, since the 
noise does not affect the populations associated with the 
computational basis, the classical capacity of the phase 
damping channel (98) C{Vg) = 1. 

On the other hand the quantum capacity of a phase 
damping channel (98) is equal to 



Q{Vg) = 1 - H{p,) = 1 - H{l/2 + g/2) 



(99) 



where H is the binary entropy. We hence have, for g > 
9i QiVg) > QiVg), i.e. enhanced quantum capacity by 
memory effects. 

Markovian correlated noise can be derived from the 



mapping (91) by properly choosing the transformation 
£r ■ Consider the case in which for sufficiently big r the 
map Sr describes a decoherent process of LE where, given 
{K)l_e} an orthonormal basis of "Hlb, one has 



SrmLE{i'\)=Se,i, \Mr))LE{Mr)\ 



(100) 
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with the vectors {|'0^(t))l£;}£ being not necessarily or- 
thogonal. The condition ( [89| can then be satisfied by 
identifying o-q with one element of the selected basis (say 
Ko)l£;), and imposing \-)pi{T > te))le = \1^)le for all I. 
In this case the mapping (91 1 can be expressed in terms 
of the operators 

Al, = LE{il\Ui\l^)LE (101) 
At,^^ = LE{ij + l\U, + Mt,{T,))LE, (102) 

which act, respectively, on the Hilbert space T-Lci and 
T-Lcj+i for = 1, • • • , n — 1. They allow us to define the 
probability distribution 



and the conditional probabilities 



(103) 



Tr„ 



i.,,,} . (104) 



Using these quantities Eq. (91) can be finally expressed 
in compact Markovian form. 



Pt-n\t-n- 



xv4,„ • • • A,, Ai, A\^ A\^--- A\^, (105) 

w ith At, = AjJ^/pi; and Aj = A j^^J JPi^^,\i^. 

[D'Arrigo, Benenti and Falci[ |2007| have considered de- 
phasing channel along this line, i.e. with Markovian corre- 
lations. Specifically, they have taken the Kraus operators 
as A^ .^o = I, A« =1 — cTz, Vj. Moreover, the conditional 
probability have been assumed as 



(1 -A*)Wj 



(106) 



where /i S [0,1] measures the partial memory of the chan- 
nel. 

Clearly the classical capacity is C = 1 also in this case. 

In order to compute the quantum capacity, one evalu- 
ates the maximum of coherent information J(<I>*^"\ pf^"^) 
by choosing = /„/2". In this case S'($(")(/0*"^)) = 
S'(p'^"^) = n. Furthermore the entropy exchange equals 
S{W), where the density operator W has components 
W,,...,„,,;...,;^ = Tr(A,,...,„p(")4,^^^^0- Here W is diag- 
onal and 

S{W) = -^p{,,}log2P{.,} ^i?(Xi,...,X„), (107) 

where H{Xi, is the Shannon entropy of the col- 

lection of random variables Xi, ...,X„ (characterized by 
the joint proba bilities ...j^ ). For a stationary Markov 
chain, one has (Cover and Thomas 1991) 



hm -iJ(Xi,...,X„) = H{X2\X,) 



where go,z = (1 ~ A*)j'o,z + M are the conditional probabil- 
ities that the channel acts on two subsequent qubits via 
the same Pauli operator, and H^qq), H{qz) are binary 
Shannon entropies. Therefore, the quantum capacity is 
finally given by 



l-poH{qo) -~PzH{q^). 



(109) 



Interestingly enough, Q increases when considering an 
higher degree of memory. In particular, the memoryless 
dephasing channel capacity, Q = Qi = 1 — H{po), is 
recovered for /i = 0, while for fJ- = 1 (perfect memory) the 
channel is asymptotically noiseless, i.e. Q = 1 ( [Bowen 



and Mancini 2004 1 



A microscopic model for correlated dephasing has been 
introduced in |D'Arrigo, Benenti and Falci[ [2007] in terms 
of a spin-boson model, where quantum information is en- 
coded in a train of qubits and a single bosonic mode rep- 
resents the memory system. Lower bounds for the quan- 
tum capacity of a qubit memory channel with both cor- 
related dephasing and damping have been evaluated nu- 
merically starting from a microscopic spin-boson model 
with Jaynes-Cummings interaction in the presence of 



stron g dephasing noise (Benenti, D'arrigo and Falci 
20091). 



|Plenio and Virmani[ |2007[ |2008| have considered an- 
other model of dephasing memory channel for qubits. It 
can be traced back to the scenario introduced in ( Giovan- 



netti and Mancini 2005 ) and schematized in Fig. 3 where 



each individual information carrier (a qubit in this case) 
interacts with a correspondent environment particle, the 
correlations being established by the environment multi- 
particle state. Specifically they have considered the case 
where the two-particle (two-qubit) interaction is defined 
by a controUed-phase gate, the environmental particle 
being the controller qubit that determines which unitary 
transformation will be applied to the carrier. As a con- 
sequence the join state of the carriers gets transformed 
through mixtures of random sequences of identity and cr^ 
operators, each sequence being characterized by a (cor- 
related) probability which depends upon the diagonal el- 
ements of the environment initial state. 

The interesting feature of this model is that it allows to 
write down explicit formulae for the associated capacities 
for the channel in terms of properties of the many-body 
environment that share a close relationship with thermo- 
dynamical quantities. In particular the quantum capac- 
ity can be expressed in terms of the regularized diagonal 
entropy of the system environment, i.e. 



= 1 - hm 



S'(diag(pem,)) 



(110) 



PoH (go) +p,H{q,), (108) 



where diag(pem)) is the environmental state in the com- 
putational basis after eliminating all off-diagonal ele- 
ments [note that the coding argument used in order to 
arrive to Eq. (|110|) was been also independently shown by 
( |Hamada[ |^2002 )]. For the special case in which the ini- 
tial state of the environment is described by a classically 
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correlated many-body system, the last term on the rhs 
of Eq. (110 1 is exactly the entropy. Hence, the capacity 



is given by 



op J n-^co n 



log2 , 



(111) 



where Z„ is the partition function for n environment 
spins, and /3 is the associated inverse temperature. In 
other words, one can exploit results from classical statis- 
tical physics in order to compute the capacity, as shown 
by Eq. (111). If on the other hand, the state penv in- 



volve also quantum correlations, Eq. ( 110 1 reduces to the 



entropy obtained when every environment qubit is com- 
pletely dephased. Before computing it, one can notice 
that (Wolf et al. 2006) showed the existence of Hamil- 



tonians exhibiting quantum phase transitions and with 
ground states being Matrix Product States (MPS) in- 
volving only rank-1 matrices. Hence, it can be shown 
that the diagonal elements of such MPSs are equal to 
the probability p of microstates in corresponding classi- 
cal Ising chains. 

For a sake of simplicity, one does focus on a transla- 
tionally invariant MPS for a ID system of 2-level parti- 
cles, in the case of periodic boundary conditions. This 
environmental state is characterized by two matrices Qq 
and Qi and is given by the following expression lip) = 
Sii...,j„ ^^{Qii ■ ■ ■QiAl'i'i ■■■in) ■ Then, by dephasing 
each qubit, the resulting unnormalized state is 



E 



tr 



.k = l 



(112) 

It is possible to show that, if |ii . . . in) has I occurrences 
of and n — I of 1, and k boundaries between Os and Is 
blocks, then the corresponding diagonal elements of p are 

p{l,n-l,k) = a'b"-'c''/m{n), (113) 

with a (resp. b) being the eigenvalue of Qo ^ Qo (resp. 
Qi (8) Qi), c being the eigenvalue of (Qo <^ Qo)(Qi 
Qi)/{ab) and 9T(n) a normalization factor. 

Therefore, by exploiting this connection, one can easily 
compute the limit in Eq. ( |110[ ) by using well known many- 
body physics methods. Fig. [7] shows the case of the 
following Hamiltonian 



2{g' - l)alal+' - (1 + g) + {g ~ 1) 



,1+1 



(114) 

In this case, one knows that the ground state is a rank-1 
MPS which possesses a non-standard 'phase transition' 
at g = 0, where indeed some correlation functions are 
non-differentiable (thou gh continuous) an d the ground 
state energy is analytic (Wolf et al. 



2006) 



B. Continuous memory channels 



Diverging gradient 




g=o 



FIG. 7 Sketch of the capacity behaviour in the case of an envi 
ronment given by the ground state of the Hamiltonian ( 114 1 



Notice the divergent gradient near the 'phase transition', i.e. 
at 3 = ( |Plenio and Virmani] |2007[ ). 



e.g., photons transmitted in optical fibers, wave guides. 



free-space, etc. (Braunstein and van Loock 
ert and Pleniol '2003; Ferraro, Olivares and 



2005 Eis- 



aris 



Weedbrook et al.^ 2012) 



2005 



Hence, an important class 
of quantum channels are those acting on CV quan- 
tum systems (Caves and Drummond 1994). For de- 



scribing CV systems, one introduces a discrete set of 
bosonic oscillators, typically a set of normal modes of the 
electromagnetic field, described by the canonical pairs 
{xk,yk\k=i....,n, the generalized "position" and "momen- 
tum" , from which one defines the vector of phase-space 
variables z := (xi, yi, . . . , a;„, y„)* and the ladder oper- 
ators {afe,aj.} obeying the canonical commutation rela- 
tions [a^,aj^] — 5hk^ The state of the bosonic system 
is described by a density operator p, or equivalently by 
the associated characteristic function x(z) = Tr[py(z) 

where V{z,) = exp ■^Y.ki^k + t^Vk^aX - - i-Vkjak 
denotes the n-mode Weyl operator (i is the imaginary 
unit). 

A remarkable class of states is that of the so- 
called Gaussian states (Eisert and Plenio 2003[ Fer- 



raro, Olivares and Paris 2005). These are the states 
having a Gaussian characteristic function, x(-^) — 
exp(tm*z — iz*7z), where m is the vector of first mo- 
ments and 7 is the covariance matrix (CM) of the phase- 
space variables. The canonical commutation relations 
lead to the uncertainty relations, which in terms of the 
CM are expressed by the matrix inequality 7 — ii7/2 > 0, 
where 



-1 

1 



(115) 



Bosonic or continuous-variable (CV) systems are of- 
ten used to encode classical and quantum information. 



is the matrix representation of the phase-space canonical 
symplectic form. 
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Going back to quantum channels in CV systems, up to 
now attention has been mainly devoted to the study of 
Gaussian channels. Gaussian channels have the charac- 
terizing properties of transforming Gaussian states into 
Gaussian states. The action of a n-mode Gaussian chan- 
nels, in the Schrodinger picture, transforms a state with 
characteristic function into 



input states has been rigorously proven only for some 



X'(z) - x(^^"'z)/(z 



(116) 



where X*^"^ is a matrix inducing a linear trans- 
formation on the 2n-dimensional phase-space vector 
z, and the function /(z ) is Gaussian, i.e. / (z) = 

exp (id^^^'z 
The linear 



1z*F'^"^z) (Holevo and Werner 



2001 1. 

term proportional to the vector d*^' ac- 
counts for a translation (displacement) of the mean 
m, while the quadratic term proportional to the ma- 
trix F^") adds a term to the CM. Gaussian trans- 
formations which are also unitary are characterized 
by the property that X*-"^ is a symplectic matrix 

(i.e., preserving the canonical form, X*^"^ flX'-^'> = 
n), and = 0. A n-mode Gaussian channel 

is hence characterized by the triad {d^"'\ X^"\Y^^^). 
The composition of two Gaussian channels with as- 



sociated triads (d^"^ , x['''> , f/"^ ) and (d^"^ , X'^"' , Y^'"> ) 



yields (^("^dl"^ 



(n) 



'■I ~f"2 J ^^2 ^^1 '^^2 -"1 ^^2 T'-'2 )■ 

It follows that, by applying suitable Gaussian unitaries at 
the input and output of the channel, one can always re- 
duce the channel in a canonical form, in which d(") = 0, 
and the matrices y^"^ take a particular symmet- 

ric form. For the case of channels acting on one or two 
modes, the reduction to canonical forms allows the clas- 
sification of Gaussian quantum channels according to in- 



variance under unitary transformations ( 


Caruso and Gio- 


vannetti 


2006 


Caruso et al. 


2008 


Holevo 


2007 


I . If one 



takes the one-mode channel as a reference point, repre- 
senting a single use of the channel, a memory channel 
is characterized by a sequence of triads (d("\ y*^"^) 
such that either d(") ^ 0^^^ d^^^ or X^"' ^ X'^^\ 

Due to the fact that the carrier Hilbert space is infinite- 
dimensional it turns out that the capacity of CV channel 
can be infinite. It is hence meaningful to introduce a 
notion of constrained capacity. Most natural choices are 
to impose a constraint on the mean energy, or on the 
mean number of bosonic excitations. Given n uses of the 
channel, the constraint on the mean excitation number 
reads 



Tr(7) 



— ^ < iV-f- - 
2n - 2 



where TV is the maximum mean excitation number per 
input mode. 

For the case of Gaussian channels, taking also into ac- 
count that the energy and the excitation number are 
quadratic in the canonical operators, it is natural to 
conjecture that the optimization is saturated by Gaus- 
sian input states. However, the optimality of Gaussian 



cases ( 


Giovannetti, et al. 


Giedke 


2007). 



2004 Wolf, Perez-Garci'a and 



If one circumvents this problem and restricts the op- 
timization to Gaussian states and/or to states which 
are separable among different channel uses, the obtained 
quantities are lower bounds on the capacity. The re- 
sult of the optimization problem under the restriction of 
Gaussian input state is also called the Gaussian capac- 
ity. These lower bounds, together with upper bounds, 
were first considered in |Holevo and Werner) |2001[ Re- 
stricting the optimization to Gaussian input states pro- 
vides dramatic simplification in the calculation. In fact, 
a simple formula exists for the von Neumann entropy 
of Gaussian states (Holevo, et al. 1999). Given a n- 



mode Gaussian state with characteristic function xi^) — 
exp (tm*z — ^z*7z), its von Neumann entropy is a func- 
tion of the CM only: 



^ = E5K-l/2): 



(118) 



fc=i 



where 
9{x) 



{x + 1) log 




1) — xlogx for a;>0,,-,^„N 
for x<Qy^^> 



and {vk\k=i,....,n are the symplectic eigenvalues of the 
CM. For a n mode Gaussian state with CM 7, the sym- 
plectic eigenvalues are computed from the eige nvalu es of 
7 f^, where $7 is the symplectic form, Eq. ( 115 1. Ex- 

plicitly, the eigenvalues of 7®^^ii7 are {:lsiLv^YkZi,....n- 
Since Gaussian channels maps Gaussian input states into 
Gaussian output states, all the entropic functions in- 
volved in the calculation of capacities, as the Holevo func- 
tion, the coherent information, and the quantum mutual 
information, can be computed starting from this expres- 
sion for the von Neumann entropy. 

Among Gaussian memory channels, one can identify a 
subclass of channels for which the memory effects can be 
unraveled. That is, by applying suitable unitary encod- 
ing and decoding transformations, n uses of such chan- 
nels are unitary equivalent to the n single-mode channels 
used in parallel. By applying known results for the mem- 
oryless setting one may then compute the capacities of 
the memory channel or estimate lower bounds in terms 
of the Gaussian capacities. 

Such a unitary mapping from n uses of a Gaussian 
memory channel to n parallel uses of independent single- 



117) mode channels was first considered in (Cerf et al. 2005 



2006 Giovannetti and Mancini 


2005 


I, and then applied 
aacities of Gaussian 


for estimating the communication ca 
memory channels in several settings ( 


Lupo, Memarzadeh 


and Mancini[ 


2009 Lupo, Pilyavets and Mancini 


2009 



and Cerf, 2009). A formal definition of the class of mem- 



ory channels that can unraveled first appeared in (Lupo 



and Mancini 20101. 
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(1) 



(1) 



(1) 



FIG. 8 Unraveling of n uses of a memory channel. Each hori- 
zontal line indicates one bosonic mode, propagating from the 
left to the right. (j>^"^ denotes n uses of the memory chan- 
nel. _E'"' and D''"-' are pre-processing and post-processing 



Gaussian unitaries. 



a(i) 



's are one-mode Gaussian channels. 



For a given n, let us considered n uses of a Gaussian 
memory channel, denoted f^*-"-', associated with the triad 
memory channel can be unraveled 
if there exist unitary transformations E^"'\ D^"', act- 
ing on n modes, such that = <S)l=i<l}k\ 
that is, n uses of the memory channel are unitary equiv- 
alent to the tensor product of n independent, but not 
necessary identical, single-mode Gaussian channel. Re- 
calling that a n-mode Gaussian unitary transformation 
is identified by a triad (d("\ 5^"^ 0) where d*^") is a 
displacement vector and S'^"'' is a symplectic matrix, 
it follows that the memory channel can be unraveled 
if and only if there exist, for any n, symplectic ma- 



trices 4"'- S^d'^ such that 4"^X(")4") = ^Li^k'^ 



'd 



and 5{3"¥(")[S'{5"^]* = 0Li^fc^'- Since the applica 
tion of unitary transformations cannot change the ca- 
pacities of the channel, they can be equivalently com- 
puted for the unraveled channel, in which each input 
mode is transformed independently (although in general 
not-identically). This mapping is depicted in Fig. [sj If 
one is interested in the calculation of constrained capac- 
ities, then one has to take in account how the input en- 
ergy is changed by the encoding unitary E^'^\ The most 
relevant case is when the input photon-number remains 

unchanged, that is, I]fe=i "I«fc = -^'"^^ Sfc=i "I«fe'^''"'- 
This is the case when the encoding unitary is a so-called 
linear passive transformation. For optical realization, 
such transformations are implemented by a network of 
beam-splitters and phase-shifters (see e.g. (Ferraro, Oli 
2005 j ) 



vares and Paris 



Let us now review results on solvable models belonging 
to the two most relevant families of models of CV memory 
channels. 

Additive-noise channel, or channel with thermal-like 
noise. In the Schrodinger picture, the input state over n 
channel uses is transformed to 



L k 



L k 



mode. If P is factorized as function of ^1,...,^^, the 
channel is memoryless, otherwise, it is a channel with 
memory. This class of Gaussian channels is characterized 
by having X^") = 1. Without loss of generality, one put 
d(") = 0. For the case of memoryless thermal noise, the 
matrix is scalar, y^^^ A^th, with Nth > 0, and the 
Gaussian quantum and classical capacities can be com- 
puted exactly (Holevo and Werner. 20011. The memory 



channel considered in (Cerf et al. 



2005 



20061 belongs 



to this class. The latter is defined for two channel uses, 
represented by two bosonic modes, which are affected by 
correlated additive noise. A generalization of this model 
to the case of more than two channel uses was introduced 



in (Schafer, Karpov and Cerf 2009), based on a multi- 



mode Markovian extension of the two-mode correlated 
noise. The additive noise has a Markovian structure with 
associated conditional probability P(^k+i\(,k- ^ suitable 
choice of the conditional probability and of the distribu- 
tion of the initial noise variable i one gets a zero-mean 
Gaussian process describing the additive noise, with an 
associated triad of the form (0, 1, F^")). Here, one consid- 
ers th e mo dels introduced in (|Schafer, Karpov and Cerf 



2009 



whic 



and (Lupo, Memarzadeh and Mancini 
has the form of a block-matrix 



2009), for 



( Mil Mi2 . . . Mi„ \ 

M21 M22 ... M2n 



\M„i M„2 ... M„„y 



(121) 



with Mi. 



a(t>\'-i\Z^, and Z 



The parame- 



1 

-1 

ters CT > and (p G [0,1) determine the noise strength and 
correlations and the index £ = 0, 1 determines whether 
the noise is symmetric of anti-symmetric under exchange 
of the phase-space coordinates. For £ = one gets the 



model considered in ( Lupo, Memarzadeh and Mancini 
2009), and for £ = 1 that considered in (Schafer, Kar- 



pov and Cerf 2009 ) . The memory channel is unraveled 



by the application of unitary pre-processing and post- 
processing, which are realized as multiport beam-splitter. 
For the case £ = 0, the channel can be unraveled for any 
n. After unraveling, the channel appears as the tensor 
product of single-mode channel with thermal noise, for 
which the optimal input states are coherent states, which 
are separable among channel uses. On the other hand, 
for £ = 1, the memory channel cannot be unraveled for 
finite values of n. However, it can be unraveled in the 
limit n — > cx), hence allowing the computation of the 
channel capacities. The unraveled channel is the tensor 
product of single-mode channel with additive noise which 
is asymmetric in the phase-space. Due to this asym- 
metry, the optimal input states for the classical capac- 
ity are multimode squeezed states, which are entangled 
among different channel uses. Differently from the case of 
discrete- variable memory channels (Macchiavello, Palma 



(120) and Virmani 2004), there is no transition in these Gaus- 



where V{S,k) is the Weyl operator for the k-th input sian models: the optimal input states are either separable 
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or entangled according to the model symmetries ( Cerf et 
~d\ |2005[ |2006[ |Lupo and Mancini[ |2010[). As entangled 



states cannot be prepared locally, it is crucial to identify 
suboptimal input states that can be prepared efficiently. 



This issue was considered in (Schafer, Karpov and Cerf 



2011), where it was shown that encoding classical infor- 



mation via Gaussian matrix-product states (Adesso and 
Ericsson 2006; Schuch, Cirac and Wolf 2008), which can 



be efficiently prepared, may allow to achieve a reliable 
communication rate close to the channel capacity. Fi- 
nally, an analysis of correlated additive Gaussian chan- 
nels beyond the case of Markovian correlations was pre- 
sented in (Schafer, Karpov and Cerf 2011). 



Another family of memory channels has been first in- 
troduced in fCiovannetti and Mancini 2005), belonging 



to the class of lossy bosonic channels. Upon n uses of 
the memory channel, n input modes interact mode-wise 
with n environmental modes by a beam-splitter trans- 
formation with given transmissivity and phase. This is 
a realization of the general scheme depicted in Fig. [Sj 
where each horizontal line represents an ingoing bosonic 
mode, each vertical line an environmental mode, and 
the boxes are the beam-splitters. Clearly, no correla- 
tions exist if the environmental modes are in a prod- 
uct state. In particular, a memoryless Gaussian chan- 
nel is obtained for a product of identical Gaussian states 
(Pilyavets, Lupo and Mancini 2009j), e.g. for the case of 
the vacuum state the memoryless purely l ossy bosonic 
channel is recovered (Giovannetti, et al 



2004 Wolf, 



Perez-Garcfa and Gicdkc, 2007j ), and for a thermal state 
one recovers the memoryless lossy bosonic channel with 
thermal noise (Giovannetti, et al. 2010). On the other 



hand, memory effects arise whenever there are correla- 
tions among the environmental modes. The environmen- 



tal state considered in (Giovannetti and Mancini 2005) 



is a n-mode squeezed vacuum, or a squeezed thermal 
state, in which multimode squeezing introduces correla- 
tions among channel uses, described by a density opera- 
tor of the form: 



SpthS^ , 



(122) 



where pth is a n-mode thermal state, with mean photon 
number equal to A'th, and 



S = exp 



(123) 



with {vk, vj,} being the canonical ladder operators associ- 
ated to the fc-th environment mode, n uses of the channel 
are associated to a triad (0, y^, (1 — 77)7cnv)i where 7onv 
is the CM of the n-mode environment state. The CM of 



the state in Eq. ( 122 ) is of the form 



7onv = ( A^th + ^ ] SS^ 



(124) 



where S is the symplectic matrix being the phase-space 
representation of the unitary in (123). As first noticed 




I 



'n-1 



m 



FIG. 9 (Color online) Left: 



memory channel (Lupo, Giovannetti and Mancini 



a single use of the lossy bosonic 
20T0l. 



Right: the n-fold concatenation of the memory channel: pho- 
tons entering in the fc-th input mode au can only emerge in 
the output ports hy with k' > k. 



in (Giovannetti and Mancini 2005) this channel model 



is unitary equivalent to the tensor product of indepen- 
dent, but not-identical, one-mode lossy bosonic Gaus- 
sian channel, that is, it can be unraveled. |Giovannetti| 
and Mancini[ |2005| provided upper and lower bounds on 



the classical capacity of the memory channel in terms of 
the capacity of the memoryless single-mode lossy bosonic 
channel after replacing the input mean photon-number N 
with an effective values N > N (for obtaining the upper 
bound) and N' < N (for obtaining the lower bound). It 



has been shown in ( Lupo, Pilyavets and Mancini 2009| ) 



using a class of environmental states (e.g. that intro- 
duced in (Pilyavets, Zborovskii and Mancini 2008 1), that 



there exists a range of parameters for which the memory 
channel can be unraveled with the use of unitary pre- 
processing transformation which do preserve the input 
energy. If the environment is in a pure state, a closed 
formula has been obtained for the constrained classical 
capacity: 



Cn = givN + (1 - v)M] , 



(125) 



where is a known function of the parameters of the en- 
vironmental state. The same approach has been used in 



(Lupo, Pilyavets and Mancini 2009) to evaluate bounds 



on the quantum and entanglement-assisted classical ca- 
pacity. A different kind of model of lossy bosonic channel 
with memory was introduced in (Lupo, Giovannetti and 
Mancini[ 2010). In this model the correlations among 



channel uses are not induced by the environmental corre- 
lations, but arise directly from the structure of the chan- 
nel. As depicted in Fig. |9) the action of the channel 
upon n uses is defined by the concatenation of n iden- 
tical unitary transformation coupling the input systems 
ai, 02, a„ with a collection of local environments 
ei, 62, ... e„ and the memory system m. Both the lo- 
cal environments and the memory system are represented 
by bosonic modes. Specifically the evolution of the A:-th 
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input mode is obtained by a concatenation of two beam- 
splitter transformations, the first with transmissivity e 
and the second with transmissivity 77, as shown in Fig. [9] 
This results in a non-anticipatory intersymbol interfer- 
ence channel (see Sec. III.B.l and III.B.3), in which pre- 
vious input states affect the action of the channel on the 
current input (Bowen and Mancini 2004). By varying 



tions ( Breuer and Petruccione 2002 1 , respectively based 



the transmissivity parameters, the model may be reduced 
to a memory less lossy bosonic channel (Giovannetti, et 



on the assumption of weak system-environment coupling, 
and that the correlation times of the reservoir are much 
smaller than the typical time scale of the system dynam- 
ics. Hence, a different terminology is commonly used in 
this context, where the term Markovian becomes synony- 
mous of memoryless. These approximations lead to the 
master equation 



al. 2004 ) (the input Ofc only influences the output 6fc ) , or 



to a channel with perfect memory ( Bowen and Mancini 



di 



pit) = Cp{t) 



2004) (all Ofc interacts only with the memory mode mi). 
For specific values of the parameters (that is, ?7 = 0, 
e = 1) j>(") descri bes a quantum shift channel (Bowen 



with the time-independent Liouvillian 
Gorini-Kossakowski-Sudarshan-Lindblad 



and Mancinil [2004| , where each input state is replaced by dGorini, Kossakowski and Sudarshanj |1976| [Lindblad 



the previous one. Such a model was extended in (Lupo, 



Giovannetti and Mancini 2010b ) to encompass memory 



effects in linear amplification processes. 

In general, the calculation of the (Gaussian) capaci- 
ties of these models of Gaussian memory channels can 
be divided in four steps: first the channel is mapped 
into the direct product of single-mode Gaussian chan- 
nels, then the optimization of the relevant entropic func- 
tion is performed mode-wise under constrained mean ex- 
citation number, then the distribution of the mean ex- 
citation number over the input modes is optimized, fi- 
nally the asymptotic limit is considered. The optimiza- 
tion of the distribution of the mean excitation number 
leads to a quantum water filling solution for the capacity 
of the memory channel, where the way the mean excita- 
tion number is distributed over input modes is analogous 



to the way water distributes into a vessel (Cover and 



Thomas 1991). Algorithms for the optimization were 



presented in QPilyavets, Lupo and Mancin i, 2009; Schafer, 
Karpov and Cerf| "2011). The study of models of Gaus- 



sian memory channels that can be unraveled has also 
stimulated and motivated a deep analysis of the com- 
munication capacities of the single-mode (memoryless) 



Gaussian channel ( 


Lupo, et al. 2011 Pilyavets, 


Lupo 


and Mancini 


2009 


Schafer, Karpov and Cerf 


2010). 


In particular ( 


Pilyavets, Lupo and Mancini 2009) and 


(Schafer, Karpov and Cerf| 2010p provided a character- 



ization of one-mode Gaussian channels, respectively for 
the case of lossy channels and additive noise, in terms of 
the solutions of the optimization problem of computing 
the Gaussian classical capacity as a function of all the 
parameters entering the model. 



VIII. QUANTUM CHANNELS FROM DYNAMICS WITH 
MEMORY 

The formalism of quantum channels describes states 
changes according to the physical laws of quantum me- 
chanics. Often, this states changes in open quantum 
systems are described continuously in time through the 
formalism of master equations. Within this approach, 
one usually resorts to the Born and Markov approxima- 



1976) 



(126) 

given by the 
expression 



Cp{t) = -i[H,p{t)]+Y,[v^p{t)Vi - -{V^V.^p} 



(127) 

where the operators H and Va respectively describe 
the Hamiltonian and non-Hamiltonian dynamical terms. 
The solution of ( 126 ) can be written, for given initial 
condition p{tQ), as p[t) — A{t — to)p{to), where A{t — to) 
defines a one-parameter family of dynamical maps which 
are completely positive and trace preserving, i.e. quan- 
tum channels. The dynamical map is formally given by 
A{t — to) = e^^~*^^'~. It obeys the homogeneous composi- 
tion law 



A(ti)A(i2) = A(ti +t2) , 



(128) 



for ti,i2 > 0. 

A class of Markovian master equations of this form can 
be obtained as the continuous-time limit of a concatena- 
tion of identical system-bath interactions. These mod- 



els, known as collision models (Alicki and Lendi 1987 
Scarani et al| [20021 [^rhal and DiVincenzoH2000[ |ziman| 



and Buzek 2005), are defined by the iterated unitary in- 
teractions of the system Q with n identical reservoirs 
E = (ei, . . . , e„). This cascade process, depicted in Fig. 
To} defines a quantum channel of the form 



$"(Pq) = Trs [C/Qe, • • • {/Qe„ 



(129) 



where f/gej 's indicate the identical unitary transforma- 
tions coupling the system Q with the environmental sys- 
tems. It is worth noticing that this expression shares the 
same structure of Eq. ([Tt]) after replacing M with Q and 
Q with E. This analogy is also apparent by comparing 
Fig. [To] with Fig. [T^. The unitary interaction can be ex- 
plicitly written as f/gej — exp (igHQeA.t), where Hq^ is 
the system-environment Hamiltonian, g is the interaction 
strength and At is the interaction time. The continuous- 
time limit is hence defined by taking the limits At 0, 
and n — cx) in such a way that nAt — T < 00, and 
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FIG. 10 The cascade structure of a collision model, defined 
by the concatenation of identical unitaries. Figure taken from 
(Scarani tt aZ.||2002|) 




FIG. 11 The cascade structure leading to the master equa- 
tion for correlated quantum channels discussed in (Giovan- 
netti and Palma| |2012l, described by Eq. (1291. Each row 



corresponds to a single collision model (see Fig. |10[ ) , and each 
column correspond to a memory channel (see Fig. [l^) 



g^Ai — > 7 < oo. In this limit the system evolution be- 
tween time and T is described by a Markovian mas- 
ter equation of the form (127 1. A more general frame- 



work has been recently introduced by ( Giovannetti and 



Palma 2012), where a collection of m quantum systems, 
Q = {qi, . . . Qm), is considered. The system qj repeat- 
edly interacts with a corresponding collection of n en- 



vironmental systems, Ej 



according to 



Eq. (129 1. Furthermore, memory effects are included in 
this setup by cascading the environmental systems as de- 
picted in Fig. [TT] , where an internal dynamics of the 
environmental systems is modeled by introducing an ad- 
ditional CPTP map M. In the continuous-time limit, 
the overall dynamics of the multipartite system Q is still 
described by a Markovian master equation, while non 
Markovian effects appear in the reduced dynamics of the 
individual systems qj. 



A. Non Markovian master equations 

The first and simplest generalization of the dynam- 
ical equation ( 126 1 is obtained by introducing a time- 



dependent Liouvillian C{t) admitting the representation 
(1271, but with time-dependent operators, H{t) and 
Va (t) . Hence, the time-dependent equation for the dy- 



namical map A(t, to 
d 



dt 



A{t,to)=Cit)A{t,to) , A(to,io)=id, (130) 



has formal solution 



A(i,to) = T exp / C{T)dT 



to 



(131) 



where T denotes time-ordering. Differently from the 
time-homogeneous case, the explicit dependence on time 
implies that the dynamical map A(t,to) is no more a 
function of 't — to' only. Notwithstanding, it still satisfies 
the inhomogeneous composition law 



A(t,s) • A(s,to) = A(t,io) 



(132) 



for any t > s > to. The Markovian character is hence 
prese rved by the time-dependent dynamical equation in 

(poi. 



A more general approach to the modeling of the 
dynamics of open quantum systems is based on 
the Nakajima-Zwanzig projection operator technique 



( Breuer and Petruccione 


2002 


iNakajima 


1958 


Zwanzig 


1960). According to that the dynamical equati 


on IS rep- 



resented as follows: 



dt 



Pit) 



IC{t ~ u)p{u) du , pito) = po . (133) 



Here memory effects are modeled by the introduction 
of the memory kernel operator IC{t). Hence, the rate 
of change of the state at time also depends on its his- 
tory, and the Markovian setting ( 126 1 is reobtained when 
/C(r) = 2S{t)C. 

The dynamical map A(t,io) associated to the non- 
Markovian evolution ( 133 ) is a solution of 



dt 



A{t,ta)^ / dTJC{t-T)A{T,to) , A(to,io)=id. 

(134) 

The dynamical map A{t,to) appears to be a function of 
both to and ti. However, one can notice that the dynam- 
ics of an open quantum system can be always understood 
as the reduced dynamics of an isolated one which includes 
the environment. Being the unitary dynamics of an iso- 
lated system homogeneous in time, it follows that, once 
the degree of freedom of the environment are taken into 
account, the dynamical map will be only a function of the 
difference 't — to\ that is, A{t, tp) = A{t — tp). Indeed, it 



has been proven in (Chrusciiiski and Kossakowski 2010) 
that any solution of ( |134[ ) is also a solution of the time- 
dependent Markovian equation 



-A(t-to) = £(t,io)A(t,io) , A(to,to) = id, (135) 
dt 

with a time-dependent Liouvillian defined by the loga- 
rithmic derivative of the dynamical map C(t — tp) := 
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^A(< — to) ■ A~^{t — to). However, the explicit depen- 
dence of the generator on the initial time 'to' implies that 
£ is effectively non-local in time. Finally, the formal so- 
lution of ( 135 1 reads 



A(t,to) = T exp 



C{t) dT 



(136) 



which shows that A(t, to) is indeed homogeneous in time. 
However, contrary to ( 131 ), it does not satisfy the compo- 



sition law, a fact which represents a signature for memory 
effects. 



where /(t) > 0, and f{t) dt < 1 to guarantee com- 
plete positivity, one obtains the memory kernel 

K.{t) = K{t){A-\A) , 

with A being an arbitrary quantum channel, and K(t) is 
defined through its Laplace transform. 



k(s) 



C. Markovian vs non-Markovian dynamics 



B. Legitimate memory kernels 

One of the fundamental problems in the theory of non- 
Markovian master equations is to find those conditions on 
the memory kernel /C(t) that ensure that the time evolu- 
tion map A(t, to) is completely positive and trace preserv- 
ing, i.e., a quantum channel. Contrary to the Markovian 
case, a full characterization of legitimate memory kernels 
is still missing. 

120121) a class of 



The Gorini-Kossakowski-Sudarshan-Lindblad repre- 



In (Chruscihski and Kossakowski 



memory kernels giving rise to legitimate quantum dy- 
namics (quantum channels) has been provided. The con- 
struction is based on a simple idea of normalization: 
starting from a family of completely positive maps sat- 
isfying a certain additional condition one is able to 'nor- 
malize' it in order to obtain legitimate, i.e., completely 
positive and trace preserving, quantum dynamics. 

First of all, one notices that the non-Markovian master 
equation ( 134 1 is trace preserving iff Tr[]C{t)p] = for any 
density operator p. Equivalently, this condition may be 
rewritten in the Heisenberg picture in terms of the dual 
of the memory kernel as K,* {t)I = 0, where / denotes the 
identity operator. 

Let us introduce an arbitrary family completely posi- 
tive N{t) which is differentiable and satisfies initial condi- 
tion A'o = id. Note that it may be represented as follows 



N{t) = id - 



F(r)dT, 



(137) 



where F{t) = —N{t). It is possible to prove that if 
F*{t)I > there exists a family of completely positive 
maps Q(t) satisfying 

Q*{t)I = F*{t)I . (138) 

and such that the formula 

K,{s) = [Q{s) - (id - sN{s))]N-\s) , (139) 

where ~ denotes the Laplace transform, defines a legiti- 
mate memory kernel JC{t). 

As an example, if one starts from 



N{t)={l- /(T)dT) id 



sentation (126) characterizes the infinitesimal generator 
of Markovian dynamics. Nevertheless, a unique oper- 
ative definition and quantifier of (non)Markovianity is 
still lacking, and several proposal have been put forward. 
Here we briefly review two of these proposal, respectively 
based on the notion of divisibility of a quantum channel, 
and distinguishability of quantum states. 

A one-parameter family of CPTP maps A(t, 0) is said 
to be infinitesimally divisible ( Wolf and Cirac[ 2008[ ) if it 
can be written as a 



A(t + T, 0) = A(t + T, t) A(t, 0) 



(140) 



with A(t + T,t) a CPTP ma p for any t, r > 0. (|Ri- 



vas, Huelga and Plenio 2010 ) define a map to be (time- 
inhomogeneous) Markovian when it satisfies the compo- 
sition law 



A(s,t) = A{s,u)A{u,t) 



(141) 



for any s > u > t.This composition rule is indeed a char- 
acteristic trait of Markovian evolution: it is the quantum 
counterpart to the classic Chapman-Kolgoniorov equa- 
tion. Clearly, a family of CPTP maps satisfying ( 141 ) 



defines a one-parameter family of infinitesimally divisi- 
ble quantum channels (140). A characterization of in- 



finitesimally divisible quantum channels can be obtained 
in terms of the quantity 



i(t) := lim 



|idd«) A(t + e,t)P+||i - 1 



(142) 



where denotes the projector into a maximally entan- 
gled state dimens ion d and idd is the identity ma p. It 
has been proven ( [Rivas, Huelga and Plenio 2010) that 
g{t) > if and only if the map A(t, 0) is indivisible. 

One fundamental property of CPTP maps is that they 
cannot increase the trace-distance, that is. 



i?(A(t,0)(pi),A(t,0)(pi)) <D(pi,p2) 



(143) 



for any pair of initial states pi, P2- If a family of CPTP 
maps is infinitesimally divisible, the monotonicity of the 
trace distance holds true locally, that is, 



dt 



Z?(A(t,0)(pi),A(t,0)(pi))<0. (144) 
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This idea has been used to define a criterion of non- 
Mar kovianity (Breuer, Laine and Piilo 2009). Accord- 
ing to that the dynamical map A{t, O) is said to be 
non-Markovian if there exist a value of t such that 
'^{pijP2',t) > 0, for some initial states pi, p2- Physically, 
this implies a temporal increase in the distinguishability 
of the two quantum states, due to a backflow of informa- 
tion from the surrounding environment. 

These two criteria allow one to define computable mea- 
sure of non-Markovianity. A natural quantifier derived 



from the criterion of (Rivas, Huelga and Plenio 

reads 



2010) 



A/'rhp(A) 



I 



1 ' 



(145) 



where I = g{t)dt- From the criterion of (Breuer, 



quantifier 



Laine and Piilo[ 2009 ) one defines the non-Markoviamty 



A/'blp(A) = sup 



Pl,P2 Ja>0 



(146) 



( Chruscinski, Kossakowski and Rivas 2011) have 



pointed out that the relation these two criteria resem- 
bles that between separable and PPT states in entan- 
glement theory. Indeed, any family of CPTP maps 
which is Markovian according to the first criterion is 
as well Markovian according to the second one, that is, 
A/'rhp(A) = implies A/'blp(A) = 0, while the converse 
is in general not true. 



D. A solvable model 

One of the few example of non Markovian dynamics 
that are exactly solvable for their communication capac- 
ities is a single qubit coupled to an environment of (non 
interacting) qubits giving rise to dephasing channel in- 
troduced by (Arshed, Toor and Lidar 2010). Such a 
model is defined by a spin-star system of A'' -I- 1 spin- 



1/2, where the central system plays the role of the infor- 
mation carrier, while the other N spins model the sur- 
rounding environment. The Hamiltonian of the system 

is of the form H — X]!!Lo ^"''^n' where the label n = 
refers to the central spin. Moreover, an interaction term 
is introduced between the central spin and the others: 

Hi — crgj^n^i dn.'^n- The N surrounding spins are ini- 
tially assumed to be in the thermal state associated to 
their own Hamiltonian. Since H commutes with Hj such 
an initial bath state is stationary. 

Due to the form of the Hamiltonian the associated 
quantum channel, mapping the initial state of the central 
spin to its state at time t, is a dephasing channel for any 
value of t > 0. The quantum capacity of this quantum 
channel can be hence computed for all times, yielding 



Qit) 



(1 + TTn) l0g2((l + TTn) + (1 - 7rAr)l0g2(l - TTjy) 

2 

(147) 



where the function Tr^r, besides depending on time, 
strongly depends on the parameters entering the dynam- 
ical model. In turns, the behavior of the quantum ca- 
pacity as function of time is as well strongly dependent 
on the couplings parameters and on the temperature of 
the bath. For generic values, recurrences in the capacity 
as function of time are of small amplitude and quickly 
vanish. On the other hand, for commensurable values 
of the parameters the quantum capacity is a periodic as 
function of time. This feature indicates the backflow of 
information from the environment to the central spin, a 
signature of non-Markovian dynamics. This is also re- 
lated to the increased distinguishability of states pointed 
out by the BLP criterion. 



IX. SUMMARY AND OUTLOOK 

It has been seen that accounting for memory effects in 
quantum processes (channels) is a daunting task. How- 
ever, the consideration of spatial and temporal memory 
effects is becoming increasingly pressing with the con- 
tinuing miniaturization of devices and with increasing 
communication rates. For what concern the capacities 
evaluation, among the large variety of quantum memory 
channels the class that leaves hopes for an exact evalu- 
ation is that of forgetful channels. For this reason it is 
of utmost importance to derive general criteria to decide 
whether or not a given channel is forgetful. Beyond that 
it would be extremely interesting to establish when mem- 
ory effects increase the capacity of a quantum channel. 
It is also worth noticing that the effects of correlations 
among errors are in close connection with the property 
of superadditivity of the minimum output entropy ( Hast- 



mgs 



2009 1 . The possible memory induced enhancement 



of capacity of a quantum channel, looked through the 



dynamical memory model sketched in Sec. VILA can be 



seen as due to a sort of Zeno effect ( Misra and Sudarshan 



1977). In fact, by frequently inserting information car- 



riers through the channel one prevents the environment 
to come back to its stationary state after the passage of 
each of them, thus less affecting the carriers themselves. 
Whereas in the case of quantum channels arising from 
nonMarkovian dynamics like that of Sec. |VIII.D the in- 
crement of capacity can be explained by the backflow of 
information from environment to system. 

Examples of quantum channels showing memory ef- 
fects are abundant in quantum information processing. 

An unmodulated spin chain has been proposed as a 



model for short distance quantum communication (Bose 



2003). In such a scheme, the state to be communicated 



over the channel is placed on one of the spins of the chain, 
propagates for a specific amount of time, and is then re- 
ceived at a distant spin of the chain. When viewed as 
a model for quantum communication, it is generally as- 
sumed that a reset of the spin chain occurs after each sig- 
nal, for instance by applying an external magnetic field, 
resulting in a memoryless channel. However, a continu- 
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ous operation without resetting corresponds to a quan- 



tum channel with memory ( Bayat et al. 2008 ) 



Another model of a quantum channel with memory 
is the so-called one-atom maser or micromaser (Benenti, 



D'arrigo and Falci 2009 ). In such a device, excited atoms 
interact with the photon field inside a high-quality optical 
cavity. If the photons inside the cavity have sufficiently 
long lifetime, atoms entering the cavity will feel the effect 
of the preceding atoms, introducing correlations between 
consecutive signal states. 

Furthermore, the propagation of the quantum elec- 
tromagnetic field in linear dispersive media, including 
the free-space propagation and through linear optical 
systems, can be described by a quantum channel with 
memory ([Giovannetti, et al. 2004b Lupo, et al. 2011b 



Shapiro! |2009| ), where" wave diffraction introduces the 
memory effects. 

Another source of correlated noise in the propagation 
of the electromagnetic field is due to atmospheric tur- 
bulence, whose effects on the signal propagation can be 
modeled as random changes of the channel's characteris- 
tics ( Semenov and Vogel^ 2009 2010). Moreover, the de- 



coherence induced by atmospheric turbulence introduces 



cross talks ( [Boyd et a/.[|2011[ [Tyler and Boyd||2009D , i.e. 
intersymbol interference effects (see Sec. III.B.3), when 
information is encoded in the transverse degrees of free- 
dom of the electromagnetic field, e.g. the orbital angular 
momentum. 

Memory effects in quantum cryptography, that is to 
say private capacities, have been briefly addressed in 



(Ruggeri and Mancini 2007b Vasiel et al. 2011) and 



needs further explorations. Then, an analysis of mem- 
ory effects in other channel uses-configurations, like zero 
error channel capacity, channels with feedback, chan- 
nels with unknown parameters and multiuser channels, 
should be pursued. 

Non Markovian dynamics are relevant in several physi- 
cal systems characterized by the interaction with a struc- 
tured environment. Examples are in the framework of 
solid state physics, as quantum dots in photonic cr ystals 
(iMadsen et all 120111 iVats, John and Buschl 120021) , and 



in the soft matter framework as the exciton dynamics 
in the protein environment (Caruso et al.\ |2009| |Plenio 



and Huelga, 2008.; Rcbcntrost, Chakraborty and Aspuru- 



Guzik 2009 Thorwart et al. 2009). 



Finally, notice that correlated noise effects could be 
investigated not only in the framework of the quantum 
version of the second Shannon theorem, but also in the 
framework of the quantum version of the first one. In 
particular one could consider compression of correlated 
sources, or entanglement manipulation in the presence of 
correlated noise. 
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Appendix A: Tools for characterizing quantum channels 

In order to quantify how much quantum information 
is preserved by some noisy channel $ the notion of dis- 
tance measures between quantum states needs to be in- 
troduced. Although there is a variety of such measures, 
here will focus on two widely used ones: trace distance 
and fidelity. The aim is to quantify how close two quan- 
tum states, represented by the density matrices pi and 
P2 € &{T-Lq), are in the Hilbert space T-Lq. In this con- 
text, the trace distance of two quantum states pi and p2 
is defined as: 



D{pi,P2) ■■= ^\\pi - P2II1 



(Al) 



with ||e||i := TrVete being the trace-norm., or norm- 
1 of the operator Q. It is clearly a symmetric function 
of the two states, i.e., D{pi,p2) = D{p2,pi), and van- 
ishes if and only if pi — p2. Furthermore, the trace 
distance is preserved under unitary transformation, i.e., 
D{UpiW ,Up2U^) — D{pi,p2), but contractive under 
CPTP map $: 



Dmp,),^p2))<Dip,,p2) 



(A2) 



Another important quantity, playing a key role in 
quantum communication theory, is the so-called fidelity 
F. Note that it does not define a distance in &{T-Lq) 
since, for instance, it is equal to 1 when the two states 
are the same. Mathematically, the quantum fidelity be- 
tween two states pi and p2, i.e., F{pi,p2), is given by 



Jozsa 


1994 


Uhlmann 


1976 



-F(P1,P2) 



Tr 



1/2 1/2 

Pi P2 Pi 



1 2 



(A3) 



It is always in the range [0,1], equal to 1 if and only 
ff Pi = P2, and symmetric, F(pi,p2) = F{p2,pi). It 
vanishes for density matrices with orthogonal supports, 
e.g., it is vanishing for orthogonal pure states. Further- 
more, it is invariant under the action of a unitary evolu- 
tion, F{UpiW ,Up2U^) = F{pi,p2), and monotonically 
increases under CPTP map. 



F($(pi),$(p2)) >F(P1,P2) 



(A4) 
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The entanglement fidelity (Schumacher 19961 is another 



quantity which is often used to gauge how much quantum 
information is preserved under the action of a quantum 
channel $ acting on a state p G &{Hq). To define the lat- 
ter quantity, one introduces a reference system, described 
by the Hilbert space Hr, isomorphic to the Hilbert space 
Hq. a purification of p is a pure state Itpp) € Hq Hu 
such that p = TrfllV'p) ("iApl- The entanglement fidelity is 
then defined as 



F{p,^) (^p|$(|V'p)(V^p|)|Vp 



(A5) 



Given a set of Kraus operators {A'j}j=i^...^d^ for the 
quantum channel the entanglement fidelity can be 
written as follows: 



Fe(p,$)-g|Tr(pif, 

3 = 1 



(A6) 



In particular, if p is a maximally mixed state, that is, 
p = IqI dim "Hq , the entanglement fidelity reads 



Fe(/Q/dimHQ,$) 



(dimHQ)2 



|Tr(X,)r (A7) 



Finally, the mathematical notion of norm of complete 
boundness, or cb-norm, is often used in the context of 
quantum channel theory. It is defined as: 



:=: sup II $ ® id„ 



(A8) 



where id„ denotes the identity map on ©(C"), and 
is defined by 



$ ® id„ 



sup ■ 
p 



||$®id„(p)|| 



IHI 



(A9) 



This norm is also called diamond-norm, denoted || • ||o- 

All these distance measures leads exactly to the 
same definitions of channel capacities fKretschmann and 



Werner 2004, ), as defined in the following in terms of 
fidelities and error correcting codes (see Secs.|v)|Vl]). 

Another distance quantifier is the quantum relative en- 
tropy, which is introduced in the next subsection. 



Appendix B: Entropic quantities 

In the study of quantum communication, entropic 
quantities play a fundamental role in characterizing 
quantum channels in terms of their efficiency as commu- 
nication lines (Barnum, Nielsen and Schumacher, 1998). 



Let us denote by 



S{p) :=-Tr(plogp) 



(Bl) 



the von Neumann entropy of a density operator p. Then, 
given a quantum channel 4' and its input state p, there 
are three important entropic quantities related to the pair 



(p, namely, the entropy of the input state S{p) {input 
entropy), the entropy of the output state S(^{p)) (out- 
put entropy), and the entropy of exchange S{p,^). The 
entropy exchange is defined as 



5(p,<f) :=5[(<i>®id)(|^p)(7^p|)]. 



(B2) 



that is, as the output entropy of the dilated channel 
(8) id) where the input state |V'p) is a purification of p. 
Such a quantity is independent of the choice of the purifi- 
cation. By construction the entropy of exchange can also 
be expressed as the output entropy of the complementary 
channel $ of $, i.e., S{p, $) = 5($(p)). Moreover, these 
three quantities satisfy the triangle inequality 



|5($(p))-5(p,$)| <5(p) 



(B3) 



From these three entropies one can construct several 
information quantities. In analogy with classical informa- 
tion theory, one defines the quantum mutual information 
between the reference system R (which mirrors the input 
Q) and the output system Q' (Petz [2C)08[ ), as 



S{p'r) + S{p'q) - S{p'^q) 

S{p) + Smp))-S{p,^). (B4) 



Notice that an important component of I(p, <f>) is the 
coherent information 



J(p,$) :=5($(p))-5(p,<i>) 



(B5) 



The quantum mutual information I satisfies the data- 
processing inequality: 



I{p, $2 o $i) < min{I(p, $i),I(p, $2)} , 



(B6) 



which in turn implies a similar inequality for the coherent 
information, i.e.. 



J(p,$2 o *i) < Jip,^i) , 



(B7) 



where $1 and <I>2 are two generic quantum channels. 

Finally, a distance-like entropic quantity is given by 
the quantum relative entropy. In particular, the quantum 
relative entropy of pi with respect to p2 is defined as: 



S'(pi||p2) := Tr [pi(logpi - logp2)] 



(B8) 



with pi, P2 being two generic density operators. This 
quantity is always nonnegative, equal to zero if and only 
if pi = P2. A useful feature of this entropy measure is 
its monotonicity when a quantum channel $ is applied 
to both states, i.e., 

5[<i>(pi)||<I>(p2)]<5(pi||p2), (B9) 
hence, two quantum states become less distinguishable 



after the action of a noisy quantum channel $ ( Lindblad 



1975) 
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Appendix C: Decomposition for non-anticipatory quantum 
channels 

In this section, one provides an explicit derivation 



of the decomposition ( 
analysis presented in ( 


17) 


based on a generalization of 


Beckman et al. 2001 Eggeling, 


Schlingemann and Werner 2001 |Kretschmann and 


Werner 2005| Piani et al. 


2006). 



CPTP maps which describes a non-antipatory quantum 
channel. Adopting the unitary representation ([2|, for 
each n integer one can define a unitary transformation 



W, 



(n) 



g„ i3„_i ••• (ji M coupling the first n carriers to a com- 
mon environment M which allows us to write 



Itl 



with il]\f = |f7)jv/(^^| being a pure fixed state of M in- 
dependent from the input state p^'- In general the uni- 
tary couplings ^ ... m ^^"^ have a complicated 
dependence upon n: however since the channel is non- 
anticipatory they must obey the following rule 

-1) Itl 

(n-l) _ rp„ r„(")i 



Tr 



q,^,M 



qn,qn 



where 



,qi-M 

Tr, 



Pq ' = i-tq^ [Pq ] is the reduced density ma- 
trix of Pq*^ associated with the first n — 1 carriers. 
Applying this relation to a generic pure state pg ■* of 
the form \'4')q„ \4')q„-->.--- -a-i one notices that the vec- 



tors W, 



(») 



(n-l) 



qn:qn-l:--- .91, Ml 



I,--- ,qi 

>«„ l'/'>g„-i,--- ,<?! ® I^)m and 
(8) \^)m are both pu- 



rification of the state ^^"~^^{p'q ^■*) constructed on the 
ancillary system formed by qn and M. Therefore there 
must exist a unitary transformation Uq^M acting on the 
latter which satisfies the identity 

\n)M = C/9„m<::'^..,,„mI^^>m .(C2) 



W, 



(n) 



' qn.qn-i,--- ,qi,M 
Iterating this n times yields 



(n) 



qn,qn 



,qi,M 



\^}m — Uq^^MUq^_^M ■ ■ ■ C^giAf |^^)m(-C3) 



which replaced into Eq. (CI) implies Eq. 



As a final remark, one notices that the above construc- 
tion is always possible if one assumes M to have dimen- 
sion (dimHq)": in particular this implies that one can 
take M to be a composite system formed by n subsystem 
wi , 7712 7 • • • , TOn initialized into the pure factorized state 
\^)m = |0)mi ® ■ ■ ■ ® |0)m„- Under this assumption the 



Appendix D: C*- and Quasi-Local Algebras 

C* — algebras are axiomatic generalizations of notions 
like continuous functions over compact spaces or bounded 
operators over Hilbert spaces. They are realized by sup- 
plying the algebraic structure of these elements (result- 
ing from the possibility to add and multiply them) by an 
adjoint operation and a norm. An easy introduction to 
C*— algebras is given in (Kadison and Ringrose[ 1997). 



Let A be an algebra over the field C, i.e. a vector 
space over C equipped with a distributive and associative 
product. Introduce an adjoint operation '\ : A ^ A 
which is an anti-linear map {\a-\- fia')^ = Xa^ + p.a'^ , with 
a^^ = a and (aa')^ = a'^a^ for a, a' E A, X, n E C Then, 
the algebra A over C equipped with f is called * — algebra. 
In turn, such an algebra will be termed C* — algebra if 
there is a norm 1 1 • 1 1 : A [0, oo) such that A is complete 
with respect to || • ||, \\aa'\\ < \\a\\ ■ \\a'\\ for all a,a' G A 
and ||a^a|| — ||a|p for all a G A. Furthermore, if a 
C*~ algebra A possesses an (identity) element 1 E A, 
such that la — a for all a E A, is called unital. Then 
on a unital C*— algebra A is possible to define a state 
as a linear functional 'ijj : A ^ C which is positive, i.e. 
i/j{a*a) > 0, and normalized, i.e. ip(l) ~ 1. 

By considering the algebra A = B{'H) on a (finite- 
dimensional) Hilbert space H together with the Hilbert- 
Schmidt product (a, b) = tr(a^6), a,b E A and the Riesz 
theorem, it can be easily seen that each state "0 on is 
represented by a unique density operator p E A, that is 
ip{a) — tr{pa) for all a E A. Notice that for infinite di- 
mensional systems there may be states not representable 
by density operators in this sense. 

A * — automorphism of a C*— algebra ^ is a one- 
to-one linear map T : A A with T{a^) = 
(T(a))^ and T{aa') = T{a)T{a') for ah a, a' E 
A. Any *-automorphism induces a family {Tz}zeJ. of 
*— automorphisms, with Tji+Z2 ^ 2^2i °'^Z2 and To = id_4, 
given by := T'^ , z E Z. 

A state V on ^ is said to be Tj— invariant, or 
T— invariant, ii ipoT = ^ holds true. Obviously 4'°Tz = 
ip for all z El. Moreover, the set of Z— invariant states 
is convex. Now comes the definition of ergodic state tp as 
an extreme point in the convex set of Z— invariant states. 

Quasi-local algebras are the proper mathematical tools 
to describe infinitely extended quantum lattice systems 



operator W 



(n-l) 



(see (Simon 1993) for a detailed review on the subject). 
Let us consider a finite-dimensional C*— algebra A and 
suppose to attach a copy An of A to each n E Z. Fur- 
thermore, let A C Z be a finite set and define A^ :— 
0nGA-^"- Then, A^ is called the algebra of observables 
belonging to all sites in A. Whenever A C A' C Z tensor- 
ing with identity 1a'\a in makes A^ sub-algebra 



only on the first n — 1 environmental components: con- 
sequently the Uq^M can be chosen to couple qj only with 
the corresponding j-th components of M. In this case 
the associated unitary representation ^ corresponds to 
the one of Fig. [I] b) . 



can be chosen to act non trivially of A^ with an embedding A^ 3 a a 



1a'\a e A^'. 



In this way the product of operators, as well as other 
algebraic operations, become well defined in the larger 
algebra A^ Thus, since the identity 1a'\a does not 

algebra 



change the norm, one can define the normed 

Jac Z:|A|<oo 



of local observables .4'°'^ := [j/^r-i-\\\^^ A^. Its norm- 
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completion 



:= A^"^ (Dl) 

is then called the quasi-local C* — algebra. 

The algebra A^ can be interpreted as the algebra of 
physical observables for a subsystem localized in the re- 
gion A c Z. The quasi-local algebra A^ then corresponds 
to the extended algebra of observables on the infinite one- 
dimensional lattice Z. In the quasi-local algebra each 
observable can be approximated uniformly by local ones, 

Any state ip on A^ induces a family of states 
{V'^}acZ:|A|<oo on A^ and A C A' implies that the re- 
striction of tp^ to A^ is equal to ip^. 

A shift T-.A^^A^on A^ is induced by 

9 a ~ a 1^ I-)- T(a) := 1^ a ~ a e A^~^^. 

Notice that the extension of T to the quasi-local algebra 
A^ is a *— automorphism. Moving from the action of the 
shift T, it is possible to introduce the notion of stationary 
state ijj on A^ when ip oT = tp holds true. The set of 
stationary states on A^ turns out to be convex. Then a 
state ip on A'^ is called ergodic if it is extremal on this 
set. 
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